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Executive  summary 


We  report  the  results  for  the  limited  scope  project  (MR-2412)  tilted  “Determining  Detection  and  Clas¬ 
sification  Potential  of  Munitions  using  Advanced  EMI  Sensors  in  the  Underwater  Environment.”  The 
project  focuses  on  characterizing  Electromagnetic  Induction  (EMI)  responses  in  the  underwater  set¬ 
ting  through  numerical  and  experimental  studies  with  the  objective  of  obtaining  a  better  understanding 
of  the  impact  of  underwater  environments  on  EMI  marine  systems  deployed  to  detect  and  character¬ 
ize  buried  munitions.  In  terrestrial  settings,  conduction  currents  can  be  ignored  in  most  soil  types 

(conductivity  ~  0.01  S/m).  The  measured  magnetic  fields  from  a  subsurface  metallic  object  in  the 
low-frequency  EMI  regime  can  be  modelled  as  if  the  object  were  in  free-space.  In  contrast,  marine 
environments  are  generally  highly  conductive  with  an  average  conductivity  of  seawater  of  around  4 

-  5  S/m.  For  a  survey  in  a  conductive  medium,  the  EMI  response  of  a  compact  metallic  target  can 
deviate  from  its  in-air  response. 

For  our  numerical  studies  we  developed  an  integral  equation  approach  for  a  layered  medium  that 
was  able  to  account  for  the  changing  conductivity  of  the  air,  marine  and  sea-bottom.  We  extended  the 
model  so  we  could  compute  both  the  background  and  scattered  field  response  from  a  highly  con¬ 
ducting  and  permeable  sphere  for  dipole  and  loop  transmitters  and  receivers.  We  earned  out  a  series 
of  synthetic  experiments  by  considering  various  factors  that  might  influence  EMI  signals,  including 
current  channeling  effects,  sea  depth,  the  size  of  a  loop,  lateral  offset  of  the  receiver,  host  conductivity, 
excitation  waveform  and  antenna  insulation. 

For  our  marine  measurements,  we  built  a  2  m  x  1  m  x  1  m  fiberglass  frame  and  encased  two 
receiver  cubes  in  epoxy  to  make  them  waterproof.  The  transmitter  loop  comprised  12  turns  of  wire 
arranged  in  a  2  m  by  1  m  rectangle.  A  24V  power-supply  was  used  to  provide  a  maximum  current 
of  11.4  A  using  a  25  Hz  base-frequency  with  a  50%  duty  cycle  waveform.  A  series  of  different 
measurements  campaigns  were  conducted  in  sea- water  depths  of  between  2  and  14  m  either  from  a 
boat  or  by  wading  into  water  adjacent  to  a  sandy  beach. 

The  response  from  a  metallic  body  immersed  in  a  conductive  medium  is  a  combination  of  the  eddy 
current  response  (ECR)  due  to  currents  generated  in  the  target  and  the  galvanic  coupling  of  currents 
through  the  body  (the  current  channeling  response,  CCR).  In  terrestrial  environments  only  the  ECR 
is  important.  Our  simulations  showed  that  the  CCR  from  a  highly  conducting  and  permeable  sphere 
embedded  in  the  air-sea-sediment  is  far  smaller  than  the  ECR  and  decays  much  faster  than  the  ECR 

at  a  rate  of  t~ 3.  For  the  interesting  time  range  of  0.1  ms  -  25  ms,  the  CCR  contributes  little  to  the 
target  signals  and  thus  can  be  ignored.  Meanwhile,  at  times  beyond  several  hundred  microseconds, 
the  ECR  response  approaches  the  value  for  the  same  object  embedded  in  free  space.  These  numerical 
observations  were  confirmed  by  measurements  of  an  insulated  and  non-insulated  105  mm  projectile 
at  a  range  of  different  receiver,  transmitter  and  object  offsets. 

When  considering  a  survey  close  to  seafloor,  we  found  that  the  decay  rate  of  the  background  re¬ 
sponse  is  affected  by  the  sea  depth,  or  equivalently  by  the  distance  of  the  sensor  from  the  air-sea 
interface.  Our  results  show  that  the  background  response  in  shallow  water  decay  faster  than  in  a 
deeper  water.  In  deeper  water  where  the  sensor  is  far  away  from  the  air-sea  interface,  the  correspond¬ 
ing  background  responses  asymptotically  approaches  the  response  of  a  half-space.  Our  simulations 
demonstrated  that  sea  depths  don’t  impact  the  scattered  field  response  from  a  buried  metallic  object. 
Measurements  conducted  in  water  depths  between  2  and  14  m  displayed  the  same  general  characteris¬ 
tics  as  our  simulations  with  the  response  in  shallower  water  (2  m)  falling  off  faster  than  the  measure¬ 
ments  at  14  m.  Observed  decay  rates  were  between  t-5/2  and  t~3. 

The  background  responses  measured  inside  and  outside  of  a  transmitter  loop  are  very  similar  and 
at  later  times  vary  little  with  offset  distance  from  the  transmitter  center.  These  aspects  of  the  offset 
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background  measurements  might  be  usefully  applied  to  calibrate  a  system  that  has  an  array  of  sensors. 

EMI  responses  depend  upon  other  geometry  survey  parameters  like  the  transmitter  loop  size  and 
stand-off  distance.  A  large  transmitter  loop  can  excite  a  strong  scattered  response  from  a  buried  object 
but  also  induces  a  larger  background  response.  The  latter  might  be  removed  by  subtraction  of  the 
fields  from  widely  separated  receivers.  To  test  this  type  of  background  suppression  method  we  con¬ 
sidered  two  possible  configurations  of  receivers:  vertical  and  radial  with  respect  to  the  transmitters. 
The  experiments  showed  that  the  vertical  offset  configuration  can  be  effective  at  removing  the  back¬ 
ground  signals  but  only  at  later  times.  On  the  other  hand,  the  results  with  radially  symmetric  receivers 
showed  that  the  differential  total  responses  are  identical  with  the  differential  scattered  responses  for 
all  components  across  all  time  range.  The  radial  receiving  configuration  provides  an  effective  way  to 
obtain  the  scattered  response  purely  from  free  space  models. 

The  conductive  sea-water  medium  causes  dispersion  of  electromagnetic  fields  and  can  delay  and 
distort  the  primary  and  scattered  fields.  For  the  primary  field  our  simulations  showed  significant 
distortion  effects  for  an  ideal  rectangular  waveform  particularly  at  large  offset  distances.  However, 
when  the  exponential  rise  time  of  the  transmitter  loop  was  considered,  the  distortion  effects  became 
much  less  significant  and  can  be  ignored  for  the  time -ranges  and  offset  distances  typically  encountered 
in  UXO  sensing.  These  theoretical  predictions  were  supported  by  measurements  of  the  primary  field 
distortion  conducted  with  a  flux-gate  magnetometer. 

During  the  off-time  field  dispersion  exhibits  increased  delay  with  an  increase  in  host  conductivity. 
Results  with  a  simulated  marine  TEMTADS  system  showed  that  the  scattered  responses  excited  from  a 
long  pulse  of  10  ms  are  stronger  than  those  from  a  short  pulse  of  1  ms  in  particular  for  late  times.  This 
demonstrates  one  method  to  further  enhance  target  signals,  which  is  equally  effective  and  equivalent 
in  the  terrestrial  realm.  Further  tests  with  simulated  marine  TEMTADS  data  showed  that  a  terrestrial 
dipole  polarizability  model  could  be  used  to  successfully  characterize  the  scattered  data  from  a  buried 
metallic  object  in  the  conducting  underwater  environment. 

We  also  considered  the  impact  of  insulation  of  the  receiver  and/or  transmitter  on  marine  sens¬ 
ing.  Our  preliminary  numerical  investigation  showed  that  the  insulation  had  negligible  effect  on  the 
background  fields. 

Our  last  set  of  measurements  focused  on  the  spectral  content  of  noise  and  signals  in  the  under¬ 
water  environment.  We  used  a  24-bit  analog  to  digital  converter  to  sample  a  receiver  cube  at  a  rate  of 
120  kHz.  Measurements  were  conducted  on  land  immediately  adjacent  to  marine  measurements 
collected  at  a  water  depth  of  4  m.  The  underwater  measurements  with  no  transmitter  present  were 
lower  noise  across  most  of  the  frequency  spectrum  compared  to  the  on-land  measurements.  When  the 
transmitter  was  operating  and  sufficiently  damped  by  a  resistor,  the  spectrum  and  off-time  behavior  on 
land  and  in  water  were  virtually  identical.  From  these  results  we  can  conclude  that  there  are  no 
significant  differences  between  on-land  and  in-water  sensing  from  the  perspective  of  the  frequency 
content  of  the  noise  and  signal. 

Based  on  the  results  of  this  limited  scope  project,  we  conclude  that  underwater  detection  and 
characterization  of  buried  metal  using  EMI  based  sensor  is  feasible.  While  the  conductive  sea-water 
introduces  complications  to  the  measurement  process,  there  are  numerous  practical  options  for  miti¬ 
gation  of  these  effects.  The  conductive  sea-water  can  impact  the  scattered  field  from  a  buried  metallic 
object,  but  typically  only  at  very  early  times  and  at  large  receiver  to  object  offsets.  Neither  of  these 
conditions  are  commonly  encountered  in  practice.  Thus  apart  from  the  practical  considerations  of  op¬ 
erating  underwater,  the  modelling  techniques  and  methods  that  have  been  successfully  demonstrated 
in  terrestrial  environments  can  be  utilized  for  marine  detection  and  characterization.  The  principle 
remaining  challenge  is  the  development  and  deployment  of  a  practical  and  effective  set  of  hardware 
for  marine  EMI  sensing. 
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Eddy  Current  Response 

EMI 
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m 
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Strategic  Environmental  Research  and  Development  Program 
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TEM 

Transient  Electromagnetic  Induction 
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Transient  Multi-sensor  towed  array  detection  system 

Tx 
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Unexploded  Ordnance 

uw 

UnderWater 

VED 

Vertical  Electric  Dipole 
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1  Background  and  Objective 


Hazardous  ordnance  items  are  present  along  coast-lines  and  in  rivers  and  lakes  in  waters  shallow 
enough  to  cause  concerns  for  human  recreational  and  industrial  activities.  The  presence  of  water 
makes  it  difficult  to  detect  and  remove  these  hazardous  legacies  induced  from  wars,  military  train¬ 
ing  and  deliberate  disposal  (Schwartz  and  Brandenburg,  2009).  A  number  of  techniques  have  been 
proposed  to  detect  and  characterize  Unexploded  Ordnances  (UXO)  and  discarded  military  munitions 
(DMM)  in  the  underwater  environment  including  acoustic  waves,  magnetometery,  and  electromag¬ 
netic  induction  (EMI)  (SERDP- Workshop,  2013;  Keenan  et  ah,  2010;  Shubitidze,  2011). 

In  recent  year's,  terrestrial  munitions  response  has  seen  significant  improvements  in  our  capability 
to  discriminate  MEC  items  from  benign  metallic  clutter.  These  advances  have  been  primarily  driven 
by  the  development  of  next-generation  EMI  sensors  designed  to  interrogate  small,  near-surface  targets. 
The  research  reported  in  this  project  concerns  underwater  sensing  using  EMI  which  is  distinct  from 
the  terrestrial  setting  in  several  respects  including  positioning  requirements  and  techniques,  noise  en¬ 
vironment  and  practical  constraints  on  deployment  of  sensor  systems.  For  marine  sites  an  even  more 
fundamental  difference  between  the  two  settings  is  the  conductivity  of  the  background.  In  terrestrial 
settings,  conduction  currents  can  be  ignored  in  a  typical  UXO  contaminated  soil  environment  (con¬ 
ductivity  around  ~  0.01  S/m)  (Shubitidze  et  ah,  2004).  Thus  the  measured  magnetic  fields  from  a 
subsurface  metallic  object  in  the  low-frequency  EMI  regime  can  be  modelled  as  if  the  object  were  in 
free-space.  In  contrast,  marine  environments  are  generally  highly  conductive  with  an  average  conduc¬ 
tivity  of  seawater  of  around  4  ~  5  S/m.  For  a  survey  in  a  conductive  medium,  the  EMI  response  of 
a  compact  metallic  target  can  deviate  from  its  in-air  response.  Norton  et  al.  (2005)  studied  the  elec¬ 
tromagnetic  induction  response  to  conducting  anomalies  of  spheroidal  shape  embedded  in  a  weakly 
conducting  host  medium.  Shubitidze  et  al  under  the  SERDP  MR- 1632  project  (Shubitidze,  2011;  Shu¬ 
bitidze  et  ah,  2008,  2009)  used  the  method  of  auxiliary  sources  (MAS)  to  study  the  scattered  responses 
of  a  metallic  object  placed  in  a  uniform  conductive  medium  (Shubitidze  et  ah,  2002,  2004).  Their  re¬ 
sults  showed  that  coupling  effects  between  the  metallic  objects  and  the  conductive  surrounding  can  be 
important  at  early  times  (or  at  high  frequencies)  but  that  the  scattered  responses  at  late  times  are  close 
to  those  in  free  space. 

One  limitation  of  this  previous  work  is  that  any  difference  between  the  conductivity  of  ocean 
sediments  and  sea-water  was  ignored.  Depending  upon  porosity,  sea  sediments  may  have  a  conduc¬ 
tivity  that  is  generally  less  than  1.0  S/m  (Chave  et  ah,  1991).  Given  the  presence  of  the  sea-sediment 
conductive  interface,  we  speculate  that  a  simple  uniform  conductive  space  (Norton  et  ah,  2005;  Shu¬ 
bitidze  et  ah,  2009)  may  not  provide  an  effective  model  of  the  EMI  responses  of  a  buried  object  when 
measurements  are  taken  close  to  the  seafloor.  In  addition,  we  lack  understanding  about  how  the  back¬ 
ground  responses  in  a  conducting  medium  might  mask  the  scattered  responses  of  a  target  of  interest. 

Within  the  project  we  aimed  to  establish  the  limitations  and  potential  for  detection  and  classifica¬ 
tion  using  EMI  data  in  a  variety  of  underwater  environments  through  the  use  of  modeling,  simulation 
studies  and  empirical  measurements.  Our  research  has  aspects  in  common  with  SERDP  Project  MR- 
2409  (Bell  et  ah,  2016)  where  the  principal  focus  is  on  empirical  measurements  of  the  background 
and  target  EMI  response  in  the  controlled  environment  of  a  large  salt  water  tank.  In  addition,  some  of 
the  research  questions  addressed  in  this  project  were  motivated  by  measurements  reported,  or  results 
obtained  in  (Bell  et  ah,  2016). 

In  the  first  part  of  this  report  we  describe  the  development  of  a  numerical  model  that  accounts 
for  the  influence  of  the  conductive  sea-water  medium  and  the  air-water  and  water-seafloor  interfaces 
on  EMI  sensing.  We  then  expand  the  model  so  it  can  be  used  to  compute  the  scattered  response 
from  a  compact  metallic  object  due  to  currents  induced  within  the  object  (the  Eddy  Current  Response, 
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ECR)  and  galvanically  coupled  through  the  object  (the  Current  Channeling  Response,  CCR).  Using 
the  model  developed,  we  then  report  on  a  series  of  synthetic  experiments  designed  to  help  understand 
a  number  of  factors  germane  to  underwater  EMI  sensing.  These  include:  (i)  the  impact  of  sea-water 
depth  on  the  background  response;  (ii)  the  lateral  variation  of  background  response  around  a  given 
transmitter  coil;  (iii)  the  effect  of  transmitter  size  and  stand-off  on  the  background  and  scatted  ob¬ 
ject  responses;  (iv)  the  use  of  differential  measurements  to  remove  background;  (v)  dispersion  of  the 
electromagnetic  fields  in  sea-water;  and  (vi)  inversion  of  simulated  data.  Motivated  by  some  measure¬ 
ments  inside  insulated  systems  by  (Bell  et  al.,  2016)  we  then  developed  a  model  to  understand  the 
response  of  a  Vertical  Mangetic  Dipole  (VMD)1  embedded  in  an  insulated  cavity.  In  the  final  part  of 
the  report  we  describe  a  number  of  empirical  measurements  which  were  conducted  either  from  a  boat 
or  from  the  shore,  in  water  depths  between  2  and  14  m.  We  used  a  prototype  EMI  sensor  comprising 
two  receiver  cubes  and  a  1  m  by  2  m  transmitter  loop.  These  empirical  measurements  augment  the 
theoretical  work  of  the  proceeding  sections. 

We  expect  that  our  results  will  offer  a  solid  basis  for  future  work  of  1)  adapting  and  developing  the 
data  processing,  analysis  and  interpretation  strategies  for  marine  EMI  sensor  data;  2)  designing  and 
deploying  advanced  underwater  EMI  sensor  systems;  3)  improving  the  ability  to  characterize  UW  sites 
contaminated  by  UXO  and  performing  cost-effective  detection  and  classification  in  the  UW  realm. 


1 A  VMD  is  a  simple  approximately  to  a  physical  transmitter  loop 
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2  Integral  equation  technique  in  a  Layered  Medium 


2.1  Problem  formulation 

Consider  an  underwater  environment  where  the  electromagnetic  properties  are  different  vertically.  A 
simplified  representation  comprises  a  three-layer  structure  consisting  of  the  upper  half-space  (region 
1,  air),  the  middle  layer  (region  2,  seawater)  and  the  lower  half-space  (region  3,  sediment  ground),  as 
schematically  shown  in  Fig.  1.  Taking  the  positive  2-axis  downward  in  the  right-handed  coordinate 
system  and  choosing  the  origin  of  the  2-axis  at  the  ocean  surface,  we  have  the  air-sea  interface  at 
z-{  =  0  and  the  sea-sediment  interface  22.  Then  the  thickness  of  the  seawater  depth  given  as  cisea  =  22. 

The  electromagnetic  properties  for  each  horizontally  infinite  and  homogeneous  layer  are  character¬ 
ized  by  the  conductivity,  permittivity,  and  permeability,  denoted  as  (cr*,  et,  //,).  et  =  eo er;,  p*  =  /io/Ai, 
where  i  =  1,2,3  and  eri  and  /jr,  are  the  associated  relative  permittivity  and  permeability  constants 
with  eo  =  8.854  x  10  12  (farads/meter)  and  //q  =  4tt  x  10“  '  (henries/meter).  Assuming  non-magnetic 
sediments,  the  magnetic  permeability  of  each  layer  can  be  set  to  that  of  free  space.  For  time-harmonic 
electromagnetic  fields  with  where  oj  is  an  angular  frequency,  t  time,  and  j  =  v/— T,  we  use  the 
complex  conductivity  a  defined  as  <7  =  a  +  jut  to  collectively  describe  conductivity  and  permittivity. 
Then  the  squared  wavenumber  is  given  as 

k2  =  (1) 

Under  the  quasi-static  approximation  of  low  frequencies  (Ward  and  Hohmann,  1988),  displacement 

Region  1 


Figure  1 :  Geometric  configuration  of  the  problem.  An  underwater  environment  is  represented  by  a 
three-layer  structure  that  consists  of  the  upper  half-space  (region  1,  air),  the  middle  layer  (region  2, 
seawater),  and  the  lower  halfspace  (region  3,  sediment  ground).  Loop  transmitter  and  receivers  are 
immersed  in  the  sea  to  interrogate  a  metallic  object  buried  below  seafloor. 
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currents  are  negligible  (i.e.,  we  <C  a )  and  the  conductivity  dominates  in  the  computation.  In  this  case, 
the  squared  wavenumber  in  Eq.  (1)  becomes  purely  imaginary  and  electromagnetic  fields  behave  in  a 
diffusive  manner. 

As  with  the  terrestrial  case,  we  assume  that  an  array  of  transmitters  and  receivers  can  be  deployed 
in  the  sea  to  interrogate  a  metallic  target  buried  in  underwater  sediments.  In  the  marine  setting,  the 
responses  to  be  recorded  are  not  only  generated  from  targets  of  interest  but  also  the  diffusive  inter¬ 
action  with  the  conductive  surroundings.  According  to  electromagnetic  scattering  theory  (Ward  and 
Hohmann,  1988;  Chew,  1999),  the  measured  responses  in  the  seawater  are  described  by 

H  (t)  =  Hinc(f)  +  Hsca(£),  (2) 

where  the  total  responses  Hv  =  [Hx  Hy  Hz\  are  the  sum  of  the  associated  incident  fields  Huic(f)  in 
the  conducting  layered  structure  and  the  scattered  fields  Hsca(f )  due  to  a  buried  target  upon  excitation. 
To  compute  the  Hinc(t)  and  Hsca(f),  we  use  an  integration  equation  technique  (Wait,  1970;  Chew, 
1999). 

2.2  Field  representation  for  fundamental  sources 

First,  we  describe  the  fundamental  solutions  of  electromagnetic  fields  (Green’s  tensors)  when  the 
elementary  sources  like  electric  and  magnetic  current  dipoles  are  used  to  generate  exciting  fields  in  a 
layered  medium.  The  field  expressions  are  formulated  in  the  wavenumber  and  frequency  domains  and 
can  be  numerically  transformed  into  the  corresponding  transient  responses  through  the  digital  filtering 
techniques  reported  in  Anderson  (1982). 

Consider  the  general  case  where  both  electric  and  magnetic  sources  are  present  in  a  medium.  In 
terms  of  the  magnetic  and  electric  potentials  A  and  F  (Harrington,  2001;  Ward  and  Hohmann,  1988), 
the  electric  and  magnetic  fields  can  be  expressed  as 


E  =  —V  x  F  —  jufiA  +  ^V(V  •  A) 

1  "  (3) 
H  =  V  x  A  —  <tF  + - V(V  •  F) 

JW/i 

where  V  is  a  vectorial  differential  operator:  J^x  +  +  Jpiz,  and  x,  y,  z  denote  the  three  unit 

coordinate  vectors.  Both  A  and  F  satisfy  the  inhomogeneous  Helmholtz  equation 

V2A  +  k2  A  =  -  J 

9  „  ,  (4) 

V2F  +  k2F  =  -M 

where  J  and  M  represents  electric  and  magnetic  sources,  respectively.  Hence  finding  H  and  E  fields 
for  arbitrary  EM  sources  reduces  to  solving  the  governing  equation  (4)  for  the  vector  potential  fields 
A  and  F. 

For  elementary  sources:  a  magnetic  current  dipole  and  electric  current  dipole  oriented  along  a  unit 
vector  '0  with  the  strength  prn  and  pe,  we  can  have  M(r)  =  rippTnd(r  —  r')  and  J(r)  =  if>pe5( r  —  r') 
in  Eq.  (4).  5  is  the  Dirac  delta  function,  r'  =  (x/ ,  y' .  z')  and  r  =  (x.  y.:  z)  represent  the  source  and 
field  point  coordinates,  respectively.  pm  and  pe  are  given  by  (Wait,  1970,  1982;  Ward  and  Hohmann, 
1988), 
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Pm  — 


jujfildS 

47 r 


Pe  = 


Idl 

47 r 


(5) 


where  I  is  the  current  flowing  in  a  loop  of  area  dS  or  a  linear  wire  element  dl.  When  ip  points  in 
unit  coordinate  direction  of  z  or  x  or  y,  we  say  that  these  corresponds  to  the  VMD  (vertical  magnetic 
dipole),  or  HMD  (horizontal  magnetic  dipole),  or  VED  (vertical  electric  dipole),  or  HED  (horizontal 
electric  dipole)  sources.  Applying  the  Fourier-Bessel  transform  to  field  equations  Eq.  (4)  in  a  layered 
medium  (Ward  and  Hohmann,  1988;  Wait,  1970),  we  can  obtain  a  spectral-domain  representation 
of  the  potential  fields  that  is  a  superposition  of  the  upper  and  downward  plane  waves  propagating 
in  the  ^-direction,  modulated  by  the  oscillatory  function  at  a  given  distance  between  the  source  and 
field  points.  In  the  source  region,  the  total  fields  are  written  as  the  sum  of  primary  and  secondary 
upward/downward  fields. 

When  a  VMD  is  immersed  in  the  seawater  in  Fig  1,  we  can  write  the  spectral  representation  of  the 
z-component  potential  Fjz,j  =  1,  2, 3  in  each  layer  as 


Fiz(p,z) 
F2  z(p,z) 
F’iz  (Pi  z) 


=m  /  B\e 71  21 Jo(\p)d\,  z  <  z\ 


=m 


A2e-^(*-*i)  +  B2e^z~Z2)  +  Ae-72|2-V|j  j^xp)dX,  Zl  <z<  z2, 

72  -1 


roo 

=m  /  A^e~^z~Z2J  Jo(\p)d\,  z  >  z2 

Jo 


(6) 


where  p  =  y/(x  —  x')2  +  (y  —  y')2,  7 j  =  (A2  —  &;2)1//2,  A  is  a  radial  wavenumber,  and  Jo(Ap)  is 

the  Bessel  functions  of  the  first  kind  of  order  0.  Aj  and  Bj  denote  the  downward  and  upward  field 
amplitudes  in  layer  j,  respectively.  The  expressions  in  Eq.  (6)  are  the  well-known  Sommerfeld-type 
integral.  Subscript  1,  2,  or  3  is  added  to  these  quantities  when  they  refer  specifically  to  the  upper  half 
space,  the  sea,  and  the  bottom,  respectively.  The  exponential  factors  in  the  representation  might  be 
interpreted  as  upward  and  downward  waves  that  attenuate  to  zero  at  z  =  ±00.  In  layer  2,  .4 2  and  B2 
are  the  coefficients  for  the  downward  and  upward  propagation  waves  from  the  air-ocean  and  seafloor 
interfaces,  in  addition  with  a  direct  wave  term  e~l2\z~z'\.  This  direct  term  can  be  important  for  a 
UXO  survey  where  the  transmitter  and  receiver  are  generally  close.  When  an  observation  is  close  to 
the  seafloor,  the  upward  propagation  wave  term  e72  ■2~22')  can  play  an  important  role.  Similarly,  if 
an  observation  is  near  the  ocean  surface,  the  downward  propagation  wave  term  e~72^z~2|  )  would  be 
dominant.  In  layer  1  and  3,  B\  and  A3  corresponding  with  the  upward  and  downward  propagating 
waves  in  the  air  and  the  sediment  layer.  There  can  be  no  wave  coming  from  z  =  therefore 

A\  =  B3  =  0. 

Amplitudes  Aj  and  Bj  can  be  determined  by  continuity  of  the  tangential  field  components  of  E 
and  H  at  the  interface  (Harrington,  2001;  Chew,  1999).  For  the  VMD  excitation,  the  electromagnetic 
boundary  conditions  for  F  on  layers  j  and  j  —  1  are  given  as  (Wait,  1970;  Ward  and  Hohmann,  1988) 


Fz,j— 1  —  Fzj 

1  diAj- 1  _  1  dFzJ  (7) 

Pj~i  dz  p,j  dz 


SERDP  2412  Final  Report 


5 


Nov-2016 


In  the  case  of  VED  excitation  in  the  sea  in  Fig  1 ,  the  spectral  representation  for  potential  field  A 
is  exactly  the  same  as  those  given  in  Eq.  (6)  but  amplitude  coefficients  must  satisfy  the  associated 
boundary  conditions: 


2  .Jj  1  —  -'4-ZJ 

1  dAzj_\  1  dAzj  (8) 

d j—i  dz  dj  dz 

When  a  VMD  is  in  the  lower  half  space  in  which  a  target  is  embedded,  we  can  write  the  corre¬ 
sponding  spectral  representation  of  fields  as 


r 

Fiz(p,z )  =rn  / 
Jo 


£ie7i(z— *0 j0(Ap)dA,  Z<Zl 


F2z(p,z) 

F3z(p,z) 


—m 


—m 


J  (A2e-^z~zi)  +  B2eyAz-z2)^  j0(Ap)dA,  zi  <  z  <  z2 
j™  +  ^e'73|2'Z'')  M^p)dX,  z  >  z2 


(9) 


The  same  form  of  the  spectral  representation  of  Eq.  (9)  applies  to  a  VED  source  buried  in  region  3. 
The  upper  and  downward  amplitudes  are  found  by  imposing  the  correspondingly  boundary  conditions 
of  either  Eqs.  (7)  or  (8). 

A  horizontally-directed  magnetic  (electrical)  dipole,  say  pointing  in  x  direction,  NOT  only  induces 
a  component  of  the  potential  along  the  dipole  direction,  FX(AX )  ,  but  also  a  component,  FZ(AZ),  in 
the  vertical  direction.  The  derivation  process  of  fields  is  similar  to  the  above  VMD  (VED)  cases. 
The  expressions  for  a  HMD  (HED)  excitation  are  more  involved  than  a  VMD  (VED)  case  as  we 
have  four  more  downward  and  upward  amplitudes  related  to  a  HMD  (HED)  to  be  determined  due  to 
electromagnetic  inhomogeneities.  For  the  complete  HMD-related  expressions  we  refer  the  interested 
reader  to  Song  et  al.  (2016).  Here  we  give  the  boundary  conditions  for  HMD  and  HED  that  must  be 
satisfied.  For  the  x-oriented  HMD,  continuity  of  the  tangential  field  components  of  E  and  H  at  the 
interface  leads  to  the  boundary  conditions  about  the  potential  F  (Wait,  1970;  Ward  and  Hohmann, 
1988): 

Fz,j— l  =  FZj 
1  (  dFXj-i  +  dFzj_  i  \ 

Pj- 1  \  dx  dz  J 

aj-iFxj-i  =  ajFx,j 

dFx  j— i  _  dFx  j 

dz  dz 

Similarly  for  the  x-oriented  HED,  the  boundary  cone 

Azj~i  =  Azj 

1  /  dAXtj_  i  dAZtj_i\ 

dj-i  V  dx  dz  J 

Pj— i  Axj—i  —  pjAXj 

dAXj—\  _ dAXj 

dz  dz 


1 

Pj 


dFx,j 

dx 


+ 


dF, 


Z,J 

dz 


(10) 


ditions  of  the  potential  A  hold 


^  f  dAx>j  dAzJ 

dj  \  dx  dz 


(11) 
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Figure  2:  Check  of  numerical  harmonic  fields  with  analytical  ones  for  the  observation  at  100  m  on  a 
0.01  S/m  of  homogeneous  half  space.  VMD:  (a)  Hx\  (b)  Hz.  HMD:  (c)  Hx;  (d)  Hz. 


Noticing  the  duality  between  the  boundary  conditions  of  VMD  and  VED,  HMD  and  HED  in  Eqs. 
(7)  and  (8),  Eqs.  (10)  and  (11),  we  can  easily  obtain  the  upward  and  downward  coefficients  for 
VED  and  HED  simply  by  exchanging  /ij  and  b3  in  those  corresponding  amplitude  expressions  to 
the  VMD  and  HMD.  All  the  above  improper  integrals  in  Eqs.  (6)  and  (9)  are  evaluated  through 
Anderson’s  digital  filtering  techniques  (Anderson,  1982).  Since  the  filtering  technique  was  designed 
to  have  narrow  bandlimited  filter  response  over  a  logarithmically  transformed  abscissa  range,  the 
numerical  integration  with  larger  values  of  A  is  adaptively  contained  in  the  algorithm. 

To  evaluate  the  accuracy  of  the  numerical  Hankel  transformation,  we  compute  the  harmonic  fields 
for  an  observation  at  100  m  on  a  0.01  S/m  homogeneous  half  space.  In  this  case,  analytical  expres¬ 
sions  for  the  harmonic  horizontal  and  vertical  components  of  a  VMD  and  HMD  are  given  by  Eqs. 
(4.56)  and  (4.59),  Eqs.  (4.126)  and  (4. 128)-(4. 129)  in  Ward  and  Hohmann  (1988).  Our  three-layer 
model  is  reduced  to  a  half  space  model  by  setting  Zi=  0  and  02  =  03  =  0.01  S/m.  Fig.  2  includes  a 
check  of  the  numerical  and  analytical  results  of  Hx  and  //-,  from  VMD  and  HMD  sources  which  are 
all  in  close  agreement.  With  the  same  half  space  model,  we  also  check  the  computation  of  transient 
fields  by  the  digital  filter  technique  described  below.  The  analytical  field  expressions  for  a  VMD  are 
given  by  Eqs.  (4.69a)  and  (4.70),  (4.72)  and  (4.74)  in  Ward  and  Hohmann  (1988).  Fig.  3  (a)-(b) 
shows  the  numerical  and  analytical  results  by  a  VMD  which  are  all  in  close  agreement.  For  HMD, 
there  are  no  analytical  formulae  available.  However  the  ^-component  from  HMD  must  be  reciprocal 
with  the  x-component  from  VMD.  Fig.  3  (c)  shows  the  reciprocity  by  comparing  the  numerical  results 
of  HMD  with  the  analytical  ones  of  VMD. 

Given  the  complete  integral  representations  for  the  potentials  in  the  three  layers,  the  fundamental 
field  components  in  any  layer  in  Fig.  1  can  be  obtained  through  Eq.  (3),  including  1)  the  incident 
fields  when  source  and  receiver  loops  reside  in  the  sea;  2)  the  incident  fields  that  penetrate  from  the 
source  in  the  sea  down  into  the  sediment  layer;  3)  Green’s  tensor  transmitting  a  source  in  the  ground 
upward  into  a  receiver  in  the  sea.  Next  we  present  the  formulation  to  compute  the  scattered  fields  from 
an  object  buried  in  the  sediment. 
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Figure  3:  Check  of  numerical  transient  fields  with  analytical  ones  for  the  observation  at  100  m  on  a 
0.01  S/m  of  homogeneous  half  space.  VMD:  (a)  Hx  and  dHx/dt\  (b)  Hz  and  dHz/dt.  (c)  HMD 
components  Hz  and  dHz/dt  vs. VMD  components  Hx  and  dHx/dt. 


2.3  Scattered  field  from  a  buried  metallic  object:  ECR  and  CCR 

Following  the  terrestrial  case,  we  also  consider  a  underwater  UXO  survey  where  the  dimension  of  a 
target  of  interest  is  small  relative  to  the  target-sensor  distance  such  that  the  primary  fields  around  the 
target  can  be  assumed  to  be  approximately  uniform.  Further,  for  sufficiently  low  frequencies,  it  is 
possible  to  represent  the  transient  scattering  of  a  metallic  object  by  an  equivalent  induced  dipole.  This 
concept  has  been  successfully  used  in  the  interpretation  of  land-based  EMI  sensing  of  metallic  objects 
(Das  et  ah,  1990;  Carin  et  ah,  2001;  Pasion  and  Oldenburg,  2001;  Barrow  and  Nelson,  2001;  Bell 
et  ah,  2001;  Miller  et  ah,  2001;  Pasion,  2007;  Gasperikova  et  ah,  2008;  Billings  et  ah,  2010;  Beran 
et  ah,  2013).  Norton  et  ah  (2005)  noted  that  the  scattered  field  from  a  spheroid  immersed  in  a  uniform 
conducting  medium  appear  dipole-like  if  the  target  size  is  small  compared  to  the  sensor-target  range. 
That  is,  we  use  an  induced  magnetic  dipole  to  model  the  so-called  eddy  current  response  (ECR)  due 
to  currents  generated  in  the  target.  For  a  conducting  medium,  the  perturbation  of  electric  field  due  to 
a  body  in  principle  would  produce  an  additional  EMI  response  that  is  called  the  current  channeling 
response  (CCR).  That  can  be  modelled  by  an  electric  dipole. 

For  an  isolated  object  that  can  be  treated  like  a  point-target  in  the  low-frequency  regime,  we  use 
magnetic  and  electric  polarizability  tensors  P m(uj)  and  Pe(u;)  to  characterize  the  scattering  property 
of  an  object.  Based  upon  the  Green’s  function  method,  we  express  the  scattered  field  FIsc;i  (rR;i;,  uS)  at 
*Rx  as 


Hsca(r  Rx,uj)  =  GHH(rRx,r',  w)Pm(w)H(r',  rTx,  w)  +  GHE(rRx,r' ,uj)Pe(u)F,(r' ,rTx,u), 

(12) 

where  the  product  Pm(w)H(r/,  rRx,  to)  is  the  induced  magnetic  dipole  moment  (Das  et  ah,  1990; 
Carin  et  ah,  2001 ;  Bell  et  ah,  2001)  in  a  target  by  the  illuminating  field  H(r',  ttx,uj)  from  a  transmit¬ 
ting  antenna  at  rRx-  Correspondingly,  the  product  Pe(w)E(r',  r Tx,u)  is  the  induced  electric  dipole 
moment.  GRR(rRx,  r,  cu)  specifies  the  magnetic  Green  tensor  in  the  background  relating  the  field  at 
rRx  due  to  the  target  at  r  radiating  into  the  sea;  r,  oj)  a  mixed  magnetic  Green  tensor  in 

the  background  produced  by  the  electric  sources.  The  primary  H  and  E  fields  and  Green’s  tensors 
G hh  and  G he  are  calculated  using  the  three-layered  formulae  presented  in  Section  3.1  for  vertical 
and  horizontal  electrical  and  magnetic  dipoles.  Eq.  (12)  is  a  simplified  form  of  a  full  volume  integral 
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equation  (Wannamaker  et  al.,  1984;  Ward  and  Hohmann,  1988;  Chew,  1999;  Song  and  Liu,  2005; 
Farquharson  et  al.,  2006)  where  we  replace  the  total  fields  internal  to  a  body  by  the  induced  dipole 
moments  at  a  point  target. 

For  the  above  approximation  to  be  used  in  a  underwater  setting,  we  neglect  the  possible  interac¬ 
tions  between  a  buried  object  and  the  seafloor  assuming  that  the  object  is  not  close  to  the  interface. 
This  makes  (12)  easily  implemented  to  compute  the  scattered  field  response  of  an  object  when  its 
magnetic  and  electric  polarizabilities  are  known.  For  a  general  case,  one  might  resort  to  solving  a  full 
domain  integral  equation  (Wannamaker  et  al.,  1984;  Ward  and  Hohmann,  1988;  Chew,  1999;  Song 
and  Liu,  2005;  Farquharson  et  al.,  2006),  which  can  be  challenging  for  very  high  conductivity  and 
permeablilty. 

For  a  highly  conducting  and  permeable  object,  its  conductivity  and  relative  permeability  can  be 
6  ~  7  and  2  ~  3  orders  of  magnitude  larger  than  those  of  its  surrounding  medium  such  as  the  sea- 
sediment  layers  of  interest.  Therefore,  it  might  be  reasonable  to  assume  a  relatively  weakly  conduct¬ 
ing,  non-magnetic  background  in  the  simulation.  With  this  assumption,  we  use  a  sphere  as  our  object 
of  interest  for  the  analytic  form  of  its  magnetic  polarizability  under  the  quasi-static  regime  (i.e.,  dis¬ 
placement  current  is  ignored)  (Wait,  1951;  Ward,  1959;  Grant  and  West,  1965).  Given  a  sphere  with 
a  radius  of  o,  relative  permeability  /rr,  fi  =  qo/U,  and  conductivity  er,  the  polarizabiltiy  tensor  can  be 
expressed  as 

Pm  =  — 27ra3i?mI,  (13) 

where  I  is  a  3  x  3  identity  matrix  and  the  complex  response  function  Rm  is  given  by 

^  2/rr(tanha  —  a)  +  (tanha  —  a  +  a2  tanhct) 

/ir(tanh  a  —  a)  —  (tanh  a  —  a  +  a2  tanh  a) 


and 


a  =  (iuiUcrY^a. 


Following  the  procedure  in  the  magnetic  case  where  an  alternating  uniform  magnetic  field  is  applied 
(Grant  and  West,  1965),  we  can  derive  the  electric  polarizabiltiy  tensor  of  a  conducting  and  permeable 
sphere  assuming  that  an  alternating  uniform  electric  field  is  applied,  i.e., 


Pe  =  —  27ra3(7ii?eI, 


(15) 


where  Re  the  electric  complex  response  function  is  given  by 

2&2  (tanh  a  —  a)  +  a \  (tanh  a  —  a  +  a2  tanh  a) 
<72  (tanh  a  —  a)  —  <74  (tanh  a  —  a  +  a2  tanh  a) 


(16) 


Using  the  small  argument  approximation  of  the  hyperbolic  function  tanh(-),  for  a  small  a  we  have 


2<7i  +  <7  2 


and  thus 


Pe 


— 47ra3<7i 


04  —  <72  j 
2<7l  +  <72 


Further  when  02  04,  Eq.  (17)  reduces  to 


Pe  =  47ra3<7il. 


(17) 


(18) 
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This  result  is  the  same  as  the  one  derived  for  an  electrical  dipole  moment  when  an  uniform  plane 
wave  illuminating  a  perfectly  conducting  sphere  in  Harrington  (2001);  Norton  et  al.  (2005).  For  an 
insulating  object  in  a  conducting  medium,  <72  <C  ci,  we  see  either  from  (16)  or  (17), 

Pe  =  —  27ra3(TiI.  (19) 

Thus,  the  response  to  a  spherical  void  is  half  the  magnitude  and  opposite  in  sign  to  the  response  to  a 
perfectly  conducting  sphere  (Norton  et  al.,  2005). 

Having  obtained  the  steady-state  solutions  to  the  background  and  secondary  response,  we  now 
proceed  with  the  computation  of  the  associated  transient  scattering  responses.  In  time-domain  EMI 
sensing  of  UXO,  the  data  are  in  general  collected  as  a  series  of  voltages  after  the  transmitting  current 
is  shut  off.  The  voltage  measurements  are  related  to  the  time  derivatives  of  magnetic  field  (Harrington, 
2001;  Ward  and  Hohmann,  1988).  For  an  idealized  step-off  current  excitation,  the  time  derivatives  of 
magnetic  fields  are  the  impulse  responses  measured  in  a  receiver  coil.  Denote  an  EM  field  quantity 
H (uj)  in  the  frequency  domain,  representing  either  the  background  fields  or  secondary  ones  in  our 
case.  The  corresponding  impulse  response  h(t)  with  an  implied  causality  constraint,  i.e.,  h(t)  =  0, 
t  <  0,  can  be  expressed  as  an  inverse  sine  transform  (Newman  et  ah,  1986;  Arfken  and  Weber,  1995), 

2  r°° 

h{t)  = - /  Im[T7(u;)]  sine vtdu,  (20) 

n  Jo 

where  Im[-]  is  an  imaginary  operation  over  a  complex  field  H{uj).  The  sine  transform  of  (20)  is  simi¬ 
larly  evaluated  with  the  digital  filtering  techniques  (Anderson,  1982).  For  a  general  current  excitation, 
the  transient  response  can  be  calculated  via  a  convolution  of  the  impulse  response  and  the  transmitter 
waveform  (Wait,  1982;  Song  et  ah,  2008). 

We  are  interested  to  know  whether  the  CCR  is  small  enough  to  be  ignored  compared  to  the  ECR. 
In  the  following  numerical  studies,  we  compare  the  ECR  and  CCR.  We  model  a  marine  environment 
by  using  the  electromagnetic  parameters:  o\  =  0  S/m,  fir \  =  1;  02  =  4  S/m,  //r2  =  1;  <73  =  0.4 
S/m,  fj,r 3  =  1,  which  represent  the  air,  sea,  and  the  sediment,  respectively.  The  air-sea  and  sea- 
sediment  interfaces  are  at  z\  =  0  m,  22  =  20  m.  The  seawater  depth  d-2  =  20  m.  A  steel  sphere  with 
a  radius  of  a  =  0.3  m,  a  =  106  S/m,  and  /ir  =  300  is  buried  below  the  sediment  at  a  depth  of  3.25 
m.  The  sphere  is  located  at  r  =  (0.05, 0.05,  23.25)  m  in  the  coordinate  system.  This  simulates  a  deep 
object  scenario  in  the  underwater  setting.  Transmitters  and  receivers  are  positioned  0.1  m  above  the 
seafloor.  In  the  experiment,  we  assumed  two  Tx-Rx  offsets:  (1)  a  transmitting  source  located  at  r tx  = 
(0.10,  —0.10, 19.9)  m  and  a  receiver  at  vrx  =  (—0.10,  0.10, 19.9)  m;  to  provide  a  distance  between 
the  Tx  and  Rx  points  of  0.28  m;  and  (2)  ytx  =  (1-0,  —1.0, 19.9)  m  and  yrx  =  (—1.0, 1.0, 19.9)  m; 
the  distance  between  the  Tx  and  Rx  points  of  about  2.83  m.  Figs.  4  (a)-(b)  and  (c)-(d)  show  the  ECR 
and  CCR  for  x-  and  2-  components  when  Tx-Rx  offsets  are  at  0.28  m  and  2.8  m.  We  see  that  CCR 
can  be  larger  than  ECR  earlier  than  0. 1  ms  and  decays  fast,  at  a  rate  of  t 

Our  numerical  results  are  consistent  with  the  study  (Norton  et  al.,  2005)  and  the  recent  experi¬ 
mental  results  (Bell  et  al.,  2016).  Norton  et  al.  (2005)  argued  that  the  CCR  might  be  significant  when 
sensor-target  distance  is  larger  than  the  skin  depth  in  a  host  medium  or  equivalently  when  measure¬ 
ments  are  concerned  with  earlier  times.  Given  that  the  earliest  time  interested  in  the  UXO  survey  is 
about  0.1  ms  and  with  the  frequency  of  10  kHz  and  the  conductivity  0.4  S/m  of  sea-sediment,  we 
estimate  a  minimum  skin  depth  of  8  m.  This  is  sufficiently  greater  than  the  deepest  possible  detection 
distance  2  ~  3  m  in  the  underwater.  Bearing  in  mind  these  survey  parameters  and  the  numerical  com¬ 
parison  of  ECR  and  CCR,  we  consider  that  the  CCR  or  electric  dipole  contribution  to  the  EMI  response 
is  negligible  in  our  case  and  thus  ignore  the  second  term  of  Eq.  (12)  in  the  following  simulations. 
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(a)  (b)  (c)  (d) 


Figure  4:  Comparing  eddy  current  response  and  current  channeling  response.  Tx-Rx  offset  =  0.28  m: 
(a)  x-component;  (b)  z-component.  Tx-Rx  offset  =  2.8  m:  (c)  x-component;  (d)z-component. 


2.4  Computing  fields  for  a  loop  transmitter  and  receiver 

In  EMI  sensing  of  UXO,  either  rectangular  or  circular  loop  transmitters  and  receivers  are  commonly 
used.  To  compute  the  fields  generated  from  a  transmitter  loop  immersed  in  the  sea,  we  might  simply 
do  a  surface  integration  over  a  transmitter  loop  Tx,  i.e., 

H(r,r Tx,u)  =  J  J  Hm(r,r Tx,u)dSTx,  (21) 

where  cIStx  is  a  surface  element  and  is  equal  to  dxdy,  dxdz,  and  dydz  for  a  z-,  y-,  and  x-directed 
rectangular  loop,  respectively.  Hm  denotes  the  field  at  r  generated  by  a  associated  unit  magnetic 
dipole.  In  other  words,  we  view  a  loop  as  consisting  of  many  elementary  magnetic  dipoles  and  nu¬ 
merically  convert  Eq.  (21)  into  a  discrete  summation  over  each  individual  dipole  source  contribution. 
An  alternative  way  is  to  use  the  electric-type  source  formulation  in  terms  of  the  A-potential  in  Sec. 
2.2.  The  field  of  a  rectangular  loop  in  the  sea  can  be  expressed  as  a  line  integral 

H(r,  ytx-,^)  =  /  He(r,  rTx,  u))dlTx,  (22) 

JTx 

where  He  is  a  magnetic  field  produced  by  an  electric  dipole  and  dlr-x  the  differential  element  along  the 
loop  edge.  For  a  z-directed  Tx  loop,  we  need  x-and  y-directed  electrical  dipole  elements  to  compute 
the  primary  field  H  and  E  at  a  specified  location  r.  For  a  x-directed  Tx  loop,  we  need  y-and  z-directed 
electrical  dipole  line  integrations. 

To  verify  the  loop  formulations,  we  assume  a  square  Tx  loop  of  88.62  m  x  88.62  m  sitting  on 
the  surface  of  a  half  space  of  0.01  S/m.  For  a  unit  step-off  excitation,  Fig.  5  (a)  shows  the  transient 
dHz/dt  fields  at  the  center  of  the  loop  computed  with  Eqs.  (21)  and  (22).  The  numerical  results  are 
obtained  by  using  the  midpoint  integration  rule  (100  intervals)  and  digital  linear  filter  technique.  The 
analytical  field  is  calculated  by  using  Eq.  (4.97)  in  Ward  and  Hohmann  (1988)  with  a  circular  loop 
of  radius  50  m.  The  circular  loop  has  the  same  surface  area  as  the  square  loop.  Both  analytical  and 
numerical  results  agree  well  with  each  other  with  an  error  of  less  than  2%  after  0.01  ms  (Fig.  5(b)). 
To  save  computational  time,  it  is  preferred  to  use  Eq.  (22)  performing  a  one-dimensional  integration 
along  loop  edges. 
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To  compute  the  primary  magnetic  flux  received  over  a  receiver  loop,  we  need  to  carry  out  an 
integral  over  the  surface  of  an  Rx  loop: 


$(riix,rTx,w) 


H{rRx,rTx,u)  ■  dSRx 


(23) 


where  dSRx  is  the  directed  surface  element  of  a  Rx  loop  and  H  is  given  either  by  Eq.  (21)  or  (22). 
Eq.  (23)  involves  a  quadruple  or  triple  numerical  integration.  However  it  can  be  simplified  by  using 
the  relationship  He  =  V  x  A  and  the  Stoke  theorem. 


${rRx,rTx,u)  =  H(rRx,rTx,u)  ■  dSRx 
J  Rx 


Tx  L  J  Rx 


A(rRx,rTx,u)  ■  d\Rx 


dlRx 


(24) 


where  d\Rx  is  a  directed  line  element.  The  calculation  of  the  primary  flux  passing  through  a  receiver 
loop  becomes  a  double  line  integration  along  Tx  and  Rx  loop  with  the  potential  A.  Similarly,  we 
assume  a  square  Tx  loop  of  88.62  m  x  88.62  m  sitting  on  the  surface  of  a  half  space  of  0.01  S/m.  For 
a  unit  step-off  excitation,  Fig.  5  (c)  shows  the  computed  voltages  for  a  centered  receiver  loop  of  0.1  m 
x  0. 1  m.  This  verifies  the  ,4-potcntial  based  integral  in  Eq.  (24). 

We  now  turn  to  the  problem  of  computing  the  secondary  magnetic  flux  or  voltage  induced  by  an 
object.  Considering  only  the  ECR  term  in  Eq.  (12)  and  integrating  over  the  loop  surface  oriented  in  a 
unit  direction  n,  we  have 


Hsca(rKx,w)  •  n dSRx 


n  •  G23(rRx,r',uj)dSRx 


Pm(r/,  cu)H(r/,  rTx,  w), 


(25) 


where  dSRx  is  the  area  element  of  a  receiver  loop.  Here  G23(rRx,  r',  oj)  is  specifically  used  to  rep¬ 
resent  a  magnetic  Green’s  tensor  for  a  case  where  a  source  is  in  region  3  and  a  field  point  in  region  2 
(Fig.  1). 

Recall  the  Lorentz  reciprocity  theorem  (Harrington,  2001), 


G23(rRx,r',u)  =  G32(r',rRx,u ;), 

which  states  that  the  fields  in  medium  2  for  source  in  medium  3  can  be  obtained  as  the  fields  G32 
received  in  medium  3  as  the  source  is  in  medium  2.  Thus  we  can  transform  the  term  in  the  bracket  of 
Eq.  (25)  to 


/  G23(rRx,r' ,u)ndSRx  =  /  G32(r/,  rRx,  u)ndSRx  =  /  Hm.n(r',r Rx,u)dSRx,  (26) 

J  Rx  J  Rx  J  Rx 


where  Hm>n  =  G32(r',  rRx,  ccjn  is  the  vectorial  field  generated  by  a  magnetic  dipole  source  oriented 
in  n,  i.e.,  parallel  to  the  normal  of  the  Rx  loop.  Similar  to  Eq.  (21)  that  is  for  a  Tx  loop,  we  perform 
a  surface  integration  over  a  tri-axial  Rx  loop  with  its  normal  direction  given  in  three  coordinates 
n  =  x,  y,  z.  Again,  the  2-dimensional  surface  computation  in  Eq.  (26)  can  be  converted  into  a  line 
integral  along  a  Rx  loop,  i.e., 

H(r  ,tRxilj')  —  /  HTO)n(r  ,TRx,tu')dSRx  —  /  He  n(r  ,vRx,uj')dlRx,  (27) 

J  Rx  J  Rx 
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where  He  n(r',  r^,  oj)  is  the  field  that  would  be  produced  by  the  electric  current  elements  associated 
with  a  n-directed  loop.  For  a  z-directed  rectangular  Rx  loop,  we  perform  numerical  integrations  over 
four  line  segments  along  x  and  y  directions;  for  a  x-di reeled  rectangular  Rx  loop,  we  need  y-and 
z-directed  electrical  dipole  line  integrations.  H(r',  rRx.  ui)  in  Eq.  (27)  can  be  viewed  as  the  magnetic 
field  at  the  target  location  r'  produced  by  a  n-directed  receiver  loop  as  if  a  unit  current  is  flowing. 
With  Eq.  (27),  the  secondary  flux  of  Eq.  (25)  can  be  re-written  compactly  as 


$sca(r  Rx,v)  =  H(r,  rRx ,  w)  •  Pm(r',  w)H(r',  rTx,  w),  (28) 

where  ®sca,(rRx,  u)  =  fRxHsca(rRx,uj)  ■  n  dSRx.  Note  that  the  polarizability  tensor  Pm(r',u;)  in 
Eq.  (28)  is  spatially  invariant.  The  argument  r'  in  it  is  used  to  merely  indicate  where  an  associated 
object  is  present. 

We  have  had  a  complete  description  of  the  integral  equation  technique  in  a  layered  medium  for 
loop  transmitters  and  receivers.  Next  we  present  various  synthetic  experiments  that  aim  to  improve 
our  understanding  of  the  transient  responses  that  occur  in  the  marine  UXO  sensing. 
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(a) 


x  Field  H-based  Surface  Integral 
O  Potential  A-based  Line  Integral 

io2 


10-3  IO"2  IO'1  10°  101 


Time  (ms) 
(C) 


Figure  5:  A  check  of  computing  fields  at  the  center  of  a  circular  loop  of  radius  50  m  on  a  0.01  S/m 
of  homogeneous  half  space,  (a)  Surface  and  line  numerical  integrations  and  analytical  responses, 
(b)  Calculating  primary  voltage  with  Eq.  (24).  (c)  Relative  errors  between  numerical  and  analytical 
results  in  (a). 
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3  Synthetic  Experiments 

In  the  following  numerical  studies,  we  model  the  marine  environment  by  using  the  electromagnetic 
parameters:  o\  =  0  S/m,  /vi  =  1;  02  =  4  S/m,  \iT 2  =  1;  <73  =  0.4  S/m,  \iT 3  =  1,  which 
represent  the  air,  sea,  and  the  sediment,  respectively.  The  air-sea  interface  is  set  at  z\  =  0  m  and 
the  sea-sediment  interface  22  (and  hence  dsea  the  water  depth)  can  vary  from  2  m  to  20  m  in  the 
simulation.  Fig.  6  schematically  shows  the  three  underwater  sensing  setups  where  we  will  investigate 
how  background  and  scattered  responses  behave  and  vary  with  respect  to  the  conductivity  and  depth 
of  the  seawater,  the  size  of  a  transmitter  loop,  sensing  height  and  offset.  In  addition,  we  will  study  the 
field  dispersion  with  respect  to  pulse  width  and  shape  of  an  excitation  current  and  field  characteristics 
of  an  insulated  VMD.  We  also  invert  simulated  marine  TEMTADS  2x2  scattered  data. 

3.1  Effects  of  sea  depths  on  responses 

First,  we  test  how  the  background  and  scattered  responses  will  be  affected  by  increasing  sea  depth 
dgea  from  2  m  to  20  m  while  other  parameters  are  fixed.  In  Fig.  6  (a),  a  square  Tx  loop  of  0.35  m 


Air  Air 


(a) 


Air 


(c) 


Figure  6:  Schematic  setups  of  Tx  and  Rx  loops  immersed  in  the  seawater  for  synthetic  experiments, 
(a)  Studying  responses  by  varying  sea  depths  r(sea  =  2,  5, 10,  20  m.  (b)  Studying  responses  by  varying 
xrx  =  0.4  m,  -0.4  m,  -1.2  m,  and  -1.7  m.  dsea  =  6  m.  (c)  Studying  the  responses  by  varying  sensor 
height  h  =  0.5  m,  1.0  m,  and  1.5  m  and  using  the  Tx  loop  size:  0.5  m  x  0.5  m,  lrnx  1  m,  and  2  m 
x  2  m.  dsea  =  10  m.  The  conductivity  of  seawater  is  cr2  =  4  S/m  and  the  sediment  03  =  0.4  S/m. 
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Figure  7:  Effects  of  sea  depths  on  the  responses  with  the  setup  in  Fig.  6  (a).  Background:  (a)  dHx/dt ; 
(b)  dHz/dt.  Scattered:  (a)  dHx/dt\  (b)  dHz/dt.  A  square  Tx  loop  of  0.35  m  x  0.35  m  and  a  small 
tri-axial  cube  of  0.1m3  are  used.  The  center  of  the  Tx  is  at  r tx  =  (—0.2,  0.2,  h)  m,  the  center  of  the 
Rx  at  yrx  =  (0.2,  —0.2,  h)  m,  sensor  height  h  =  0.2  m.  The  sphere  is  with  the  radius  of  a  =  0.30  m, 
a  =  3.5  x  10' ;  //r  =  300.  Its  center  is  at  r  =  (0.05, 0.05,  dsea  +  2)  m.  The  half  space  represents  a 
sea-sediment,  approximately  equivalent  to  a  very  deep  sea  scenario. 


x  0.35  m  and  a  small  tri-axial  cube  of  O.lrn3  are  immersed  in  the  sea.  The  sensor  height  is  h  =  0.2 
m  above  the  seafloor.  The  center  of  the  Tx  is  at  ytx  =  ( — 0.2, 0.2,  dsea  —  h),  the  center  of  the 
Rx  at  r =  (0.2,  —0.2,  dsea  —  h)  m.  A  deep  steel  sphere  is  buried  at  a  depth  of  2  m  below  the 
seafloor  with  the  radius  of  a  =  0.30  m  and  a  =  3.5  x  10'  and  =  300.  The  object  is  located  at 
r  =  (0.05,  0.05,  dsea,  +  2)  m. 

For  the  background  responses  shown  in  Figs.  (7)  (a)  and  (b),  we  see  that  after  ~  0.1  ms  both 
dHx/dt  and  dHz/dt  decay  at  different  rates  for  different  sea  depths.  The  shallower  the  sea  depth, 
the  faster  the  decay  of  the  background  response.  With  increasing  sea  depth,  the  horizontal  and  ver¬ 
tical  background  responses  increase  and  asymptotically  approach  those  of  a  half-space  (gray  dashed: 
resembling  a  deep  sea  of  greater  than  50  m  water  depth)  with  decay  rate  following  t~  ~3  and  t-5/2. 
The  results  illustrate  that  the  air-sea  interface  effects  become  apparent  after  a  time  dictated  by  the 
sea  depth.  The  deeper  the  sea,  the  later  in  time  before  the  interface  effect  becomes  apparent.  For  a 
deep  sea,  the  decreasing  effect  of  the  air-ocean  interface  may  be  explained  by  the  significant  reduction 
of  the  downward  amplitude  coefficient  A2  in  Sec  2.2.  For  a  shallow  water  case,  the  presence  of  the 
air-sea  interface  diminishes  the  background  response. 

3.2  Effects  of  receiver  offset  position  on  background  response 

We  now  examine  how  the  EMI  background  and  target  signals  can  vary  with  different  Tx-Rx  offset 
and  sensor-target  distance.  In  this  test  (Fig.  6  (b)),  we  simulate  a  rectangular  Tx  loop  of  2  m  x  1 
m  immersed  in  the  sea  (dsea  =  6  m).  The  center  of  the  Tx  is  at  r tx  =  (0.0, 0.0,  dgea  —  h).  The 
measurements  are  taken  along  the  x-axis  at  xrx  =  0.4  m,  -0.4  m,  -1.2  m,  -1.7  m  and  are  marked  as 
Rxi,  Rx2,  RX3,  RX4.  The  sensor  height  is  h  =  0.5  m  above  the  seafloor. 

Figs.  8  (a)  show  the  horizontal  component  response.  Outside  the  Tx  loop,  the  response  appears 
to  shift  slightly  to  late  time  relative  to  those  inside.  The  sign  reversals  occur  twice  in  this  time  range. 
Overall,  the  background  responses  inside  and  outside  the  transmitter  loop  decay  with  a  similar  pattern. 
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Figure  8:  Background  responses  varying  with  the  Rx  offsets  in  Fig.  6  (b):  (a)  dHx/dt;  (b)  dHz/dt. 
Scattered  response  at  Rx2  in  Fig.  6  (b)  in  free  space  and  air-sea-sediment  media:  (c)  dHx/dt.  (d) 
dHz/dt.  A  z-directed  Tx  loop  of  2  m  x  1  m  and  and  a  small  tri-axial  cube  of  0.15m3  are  used.  The 
sensor  height  h  =  0.5  m.  The  center  of  the  Tx  is  at  r tx  =  (0.0,  0.0,  h)  m,  the  center  of  the  Rx  at 
?Rx  =  (®Rx5  0.0,  h )  m.  The  sphere  is  with  the  radius  of  a  =  0.05  m,  a  =  3.5  x  10' ;  p, r  =  300.  Its 
center  is  at  r  =  (0.05,  0.05,  dsea  +  0.3)  m. 


The  pattern  of  the  vertical  component  responses  in  Figs.  8  (b)  looks  different  from  the  horizontal 
ones.  The  sign  reversals  of  the  fields  only  occur  for  the  measurements  outside  the  Tx  loop  and  appear 
at  different  times  for  different  Rx  horizontal  offset.  These  can  be  explained  by  the  concept  of  the 
“smoke  ring”  that  qualitatively  describe  the  diffusion  of  the  induced  current  loop  when  a  primary 
current  is  turned  off  (Nabighian,  1979).  As  the  smoke  ring  diffuses  outward  and  its  maximum  passes 
through  an  observation  point,  the  associated  magnetic  field  flips  sign.  For  the  observation  farther 
away  from  the  Tx,  say  xrx  =  — 1.7  m,  the  corresponding  sign  changes  significantly  later  than  the  one 
at  xrx  =  — 1.2  m.  But  after  about  0.3  ms  in  this  example,  all  vertical  background  responses  either 
inside  or  outside  are  virtually  identical.  This  would  be  expected  physically  as  those  observations  sense 
the  diffusion  from  a  smoke  ring  of  ever  increasing  diameter.  These  results  indicate  that  background 
removal  for  the  vertical  field  may  be  quite  straightforward  in  the  marine  environment,  for  times  greater 
than  0.3  ms. 

For  the  same  Tx-Rx  loop  setup,  we  consider  the  scattered  responses  measured  at  Rx2  from  a 
sphere  buried  at  a  depth  of  0.3  m  below  the  seafloor.  The  sphere  is  with  the  radius  of  a  =  0.05 
m  and  o  =  3.5  x  10'  and  pr  =  300.  It  is  located  at  r  =  (0.05, 0.05,  dsea  +  0.3)  m.  Figs  8  (c) 
and  (d)  show  the  horizontal  and  vertical  fields  for  the  two  backgrounds:  free  space  and  the  air-sea- 
sediment.  The  scattered  responses  of  the  sphere  embedded  in  the  air-sea-sediment  exhibit  early-time 
sign  reversals.  After  ~  0.1  ms,  the  scattered  response  for  the  sphere  in  the  full  marine  environment 
approaches  the  response  observed  in  free  space.  These  scattered  field  observations  are  consistent  with 
the  results  for  point  transmitters  and  receivers  immersed  in  the  sea  (Song  et  ah,  2016),  other  numerical 
investigations  in  a  uniform  conducting  medium  (Shubitidze  et  ah,  2008),  and  recent  experimental 
observations  (Bell  et  ah,  2016)  in  a  large  salt  water  tank.  Namely,  the  transient  scattered  responses  of 
the  sphere  embedded  in  free-space  or  within  a  uniform  or  layered  conducting  backgrounds  all  have 
approximately  the  same  late-time  behaviour.  Again,  this  implies  that  we  can  assume  a  free-space 
background  to  compute  the  secondary  field  response. 
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3.3  Effects  of  Tx  loop  size  and  stand-off  distance 

In  this  section,  we  investigate  the  EMI  response  in  an  underwater  environment  with  a  transmitter  loop 
having  different  sizes:  0.5  m  x  0.5  m;  lmx  1  m;  and  2  m  x  1  m.  The  measurements  in  Fig.  6  (c)  are 
simulated  at  three  stand-off  distances  above  the  sea  floor:  h  =  0.5  m.  h  =  1  m,  and  h  =  1.5  m  for  each 
Tx  loop.  Figs.  9  (a)  and  (d)  shows  the  horizontal  and  vertical  background  and  scattered  responses  at 
h  =  0.5  m.  Since  the  Rx  cube  is  centered  at  the  Tx  loop,  there  is  no  horizontal  background  response. 
One  can  see  that  a  large  Tx  loop  (e.g.  2  m  x  1  m)  excites  a  stronger  scattered  response  than  with  a 
small  loop  (e.g.  0.5  m  x  0.5  m).  In  particular,  for  the  larger  stand-off  distance  of  1  m  and  1.5  m,  the 
benefit  of  increasing  scattered  response  becomes  pronounced  with  a  larger  Tx  loop  (see  Figs.  9  (b)-(e) 
and  Figs.  9  (c)-(f)).  However,  a  larger  Tx  loop  can  also  generate  a  larger  background  response  as  well. 
But  this  side  effect  can  be  mitigated  by  an  appropriate  background  removal  approach.  We  will  discuss 
this  in  the  next  section. 

When  increasing  the  sensor-target  distance,  we  know  that  the  amplitude  response  from  a  target 
will  decrease  for  a  given  loop.  There  are  other  effects  of  the  conductive  background  as  can  be  seen  in 
Figs.  9  (a)-(c)  for  the  horizontal  responses  and  Figs.  9  (d)-(f)  for  the  vertical  responses  at  the  three 
sensor  heights:  h  =  0.5  m  to  h  =  1  m  to  h  =  1.5  m.  As  the  stand-off  distance  increases  we  see  a 
“delay  propagation”  effect  in  a  conducting  medium  for  the  target  response.  This  is  most  evident  by 
looking  at  the  times  when  the  responses  change  sign.  For  the  2  m  x  1  m  loop  case,  sign  reversal  of  the 
vertical  component  occurs  at  0.004  ms,  0.009  ms,  and  0.015  ms  for  h  =  0.5  m,  h  =  1  m,  and  h  =  1.5 
m,  respectively. 

3.4  Differential  Measurements:  Background  Removal 

We  have  learned  that  the  scattered  responses  from  a  highly  conducting  and  permeable  sphere  in  a  lay¬ 
ered  conducting  medium  approach  those  in  free  space  at  late  times.  However,  practical  measurements 
are  of  total  fields  that  are  the  sum  of  the  background  and  scattered  fields  in  Eq.  (2).  To  obtain  the 
scattered  responses,  one  has  to  remove  the  background  field  from  the  observed  field.  In  terrestrial 
environments,  background  measurements  are  collected  in  target-free  regions  that  are  assumed  to  be 
identical  with  the  surroundings  of  the  target-of-interest  and  background  correction  is  a  simple  case 
of  subtraction.  In  this  section  we  experiment  with  alternative  background  mitigation  strategies  that 
utilise  differential  measurements  across  the  different  receivers. 

One  potential  configuration  in  Fig.  10  (a)  is  to  positions  4  receiver  loops  along  the  z-axis  of  a  Tx 
loop.  Sequentially,  the  4  receivers  are  marked  as  Rxl,  Rx2,  Rx3,  and  Rx4.  Their  heights  relative  to 
the  Tx  plane  are  0.0  m,  0.5  m,  1.0  m,  2  m,  respectively.  Fig.  10  (b)  shows  the  background  (dashed) 
and  scattered  (solid)  responses.  The  background  responses  measured  at  different  heights  of  Rx  are 
virtually  identical  for  times  exceeding  0.01  ms.  The  target  response  dominants  the  background  re¬ 
sponse  for  the  centered  Rx  after  ~  1  ms.  To  subtract  the  background  signals  from  the  measurements, 
we  form  the  differential  combinations  of  Rxl-Rx2,  Rxl-Rx3,  and  Rxl-Rx4.  The  corresponding  quan¬ 
tities,  e.g.,  between  Rxl-Rx4,  will  be  referred  to  as  D14ZZ-sca,  D14ZZ-bkg,  and  D14ZZ-tot  to 
represent  the  differential  scattered,  background,  and  the  total  z-eomponents,  respectively.  The  same 
notation  rule  applies  to  the  other  Rx  combinations.  Fig.  10  (c)  and  (d)  shows  the  differential  re¬ 
sponses  derived  from  the  three  combinations.  The  differential  scattered  responses  become  stronger 
than  the  differential  background  earlier  than  1  ms.  The  differential  total  response  overlaps  well  with 
the  the  differential  scattered  response  after  0.3  ms.  These  experiments  show  that  the  top-bottom  Rx 
configuration  provides  a  possible  way  to  mitigate  the  background  signals. 

Next  we  look  at  an  alternative  way  to  remove  background  signals  by  using  radially  distributed 
receivers.  As  expressed  in  the  spectral  representations  in  Sec.  2.2  the  background  EM  fields  from 
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Figure  9:  Background  (dashed)  and  scattered  (solid)  responses  varying  with  Tx  loop  size  and  sensor 
height  h  in  Fig.  6  (c).  dHx/dt:  (a)  h  =  0.5  m;  (b)  1  m;  (c)  1.5  m.  dHz/dt\  (d)  h=  0.5  m;  (e)  1  m; 
(f)  1.5  m.  The  centers  of  the  Tx  and  the  Rx  are  ytx  =  r rx  =  (0.0, 0.0,  h)  m.  The  parameters  of  the 
sphere  and  the  Rx  loop  are  the  same  as  those  used  in  Fig.  8.  Notice  that  the  geometrical  symmetry  of 
the  centered  Rx  w.r.t  the  Tx  makes  the  background  response  dHx/dt  vanished. 


both  the  dipole  and  loop  source  excitations  are  the  superposition  of  cylindrical  wave  functions  that 
propagate  vertically  in  the  form  of  e±yz  and  are  modulated  by  an  oscillatory  Bessel  function  Jo(Xp) 
or  J\  (A p)  at  a  radial  distance  p.  This  means  that,  for  a  given  excitation,  the  observed  background  fields 
would  be  the  same  at  two  different  points  as  long  as  they  are  at  the  same  height  and  keep  the  same 
distances  p  to  the  source.  Using  this  distance-dependent  EM  field  characteristic  we  can  develop  ca 
background  suppression  method.  Previously  we  used  point  transmitters  and  receivers  to  demonstrate 
background  removal  feasibility  by  arranging  two  receiver  cubes  so  they  are  radially  symmetric  with 
respect  to  a  transmitter  loop  (Song  et  ah,  2016).  Following  that  numerical  experiment,  we  conceived  a 
prototype  marine  concept  2  shown  in  Fig.  11,  which  consists  of  a  loops  with  a  size  1.8  m  x  1.8  m  and 
multiple  receivers  cubes  that  are  0. 15  m  x  0.15  m  x  0. 15  m.  The  receivers  are  arranged  symmetrically 
about  the  transmitter. 

'The  Tx-Rx  arrangement  shown  is  just  one  part  of  a  larger  multi-Tx,  multi-Rx  system 
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Figure  10:  Experiment  of  differential  measurements,  (a)  Setup  of  vertical  offsets  of  a  small  Rx  coil 
of  0.15  m  x  0.15  m  at  zrx  =  0.0  m,  0.5  m,  1.0  m,  2  m  relative  to  the  Tx  loop  of  of  2  m  x  1 
m.  (b)  Background  (dashed-lines)  and  scattered  (solid-lines)  responses  for  all  Rxs.  (c)  Differential 
background  and  scattered  responses,  (d)  Differential  total  and  scattered  responses.  The  three-layer 
parameter:  a\  =  0  S/m,  <72  =  4  S/m,  <73  =  0.4  S/m;  dsea  =  10  m.  D12ZZ-bkg,  D12ZZ-sca, 
and  D12ZZ-tot  denotes  the  differential  responses  between  the  the  background,  the  scattered,  and 
the  total  z-components  in  the  three-layered  medium.  The  same  notation  rule  applies  to  other  Rx  pairs. 
The  parameters  of  the  sphere  are  the  same  as  those  used  in  Fig.  8. 


For  illustration,  we  assume  in  Fig.  1 1  that  the  small  red  circle  in  the  figure  represents  the  projected 
position  of  a  spherical  object.  With  respect  to  this  Tx,  we  can  form  many  pairs  of  symmetric  receivers 
that  theoretically  have  the  same  EM  background  fields  due  to  the  symmetrical  distribution  of  receivers. 
For  example,  the  pairs  of  Rxl  and  Rx5,  and  of  Rx7  and  Rx9.  In  each  subplot  of  Fig.  12,  the  gray 
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Figure  1 1:  A  plane  view  of  marine  sensor  system  concept  that  consists  of  a  1.8  m  x  1.8  m  transmitter 
and  multiple  receiver  cubes  (each  0.15  m  x  0.15  m  x  0.15  m)  distributed  symmetrically  about  the 
transmitter  loop.  The  center  of  the  Tx  is  at  ytxi  =  (0.0, 0.9,  h)  m.  The  sensor  height  h  =  0.3  m.  The 
parameters  of  a  conducting  and  permeable  sphere:  a  =  0.25  m;  a  =  3.5  x  107;  //r  =  300.  Its  center 
is  at  r  =  (—0.71, 1.45,  dsea  +  2)  m  and  dsea  =  10  m.  The  small  circle  denotes  the  projected  position 
of  a  spherical  object  in  the  experiment. 

curves  represent  all  measured  responses  and  red  and  black  curves  are  used  to  highlight  the  responses 
associated  with  a  particular  Rx.  We  see  that  the  background  fields  can  significantly  mask  the  scattered 
fields  before  1  ms  or  even  later.  Yet  the  strong  background  fields  are  identical  across  Rxl  and  Rx5,  as 
well  as  across  Rx7  and  Rx9.  Conceptually,  combining  recordings  between  these  cubes  could  provide 
a  gradiometer  type  measurement  which  would  completely  remove  the  background  signals. 

Fig.  13  presents  the  differential  responses  between  the  receivers  of  Rxl  and  Rx5.  In  the  sub¬ 
plots,  D15XX-sca  denotes  the  differential  measurements  between  the  scattered  ^-components  in  the 
three-layered  medium,  D15XX-tot  the  total  x—  component  in  the  three-layered  medium.  The  same 
notation  rule  applies  to  the  y-  and  z- components  and  other  receiving  combination  like  Rx7  and  Rx9. 

Since  the  background  fields  in  the  marine  enviroment  can  be  canceled  out  across  the  two  symmet¬ 
rical  receivers,  the  three  differential  components  derived  from  the  total  fields  of  Rxl  and  Rx5  overlap 
completely  with  the  differential  ones  from  the  pure  scattered  responses  at  the  two  receivers  (Fig.  13 
(a)-(c)).  The  same  differential  results  for  Rx7  and  Rx9  can  be  achieved  for  the  three  components  (Fig. 
13  (d)-(f)).  The  experiments  show  that  the  differential  measurements  from  a  pair  of  symmetrical  re¬ 
ceivers  completely  suppress  the  background  fields  from  the  conductive  seawater.  We  know  that  within 
the  interesting  time  range  of  0.1  ms  ~  25  ms,  the  scattered  responses  in  a  layered  conducting  medium 
approach  those  in  free  space.  Thus  differential  responses  can  be  modelled  as  if  the  measurements 
occur  in  free-space.  Therefore  their  use  can  make  underwater  EMI  inversion  processing  as  simple  as 
the  land-based  case. 

3.5  Dispersion  of  Electromagnetic  fields  in  seawater 

In  the  above  modelling  and  simulation  work,  we  characterized  the  fundamental  background  and  scat¬ 
tered  responses  to  a  conducting  layered  medium  and  a  target  under  an  ideal  step-off  current  excita¬ 
tion.  In  reality,  a  realizable  current  source  generates  an  exponentially  increasing  current  flow  within 
a  transmitter  loop.  The  electromagnetic  response  due  to  a  pulse-like  excitation  has  been  known  to  be 
dispersed  in  a  conducting  medium  (Wait,  1953).  The  resulting  broadening  and  delay  due  to  conductiv¬ 
ity  and  observation  distance  may  impact  underwater  EMI  sensing  of  metallic  objects.  To  clarify  these 
questions  and  better  understand  the  propagation  of  the  transient  response  due  to  a  pulse  excitation, 
we  present  a  dispersion  analysis  with  a  VMD  or  a  loop  transmitter  embedded  in  uniform  conducting 
medium  (conductivity  of  4  S/m)  and  a  three-layered  medium  (air-seawater-sediment). 

Consider  the  transient  response  to  an  arbitrary  current  excitation  I (t)  and  a  causal  step-off  re- 
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Figure  12:  Background  and  scattered  responses  in  all  three  components  measured  from  the  marine 
system  in  Fig.  11.  Rxl:  (a)-(c).  Rx5:  (d)-(f).  Rx7:  (g)-(i).  Rx9:  (j)-(l). 


sponse  for  a  system  a(t)  (Song  et  ah,  2008),  the  field  6(f)  and  its  time-derivative  are  given  by 


t)/(t)/t 


a(0)/(f), 


(29) 


and 


db{t) 

dt 


t)I'  (r)dr  —  a(0)/,(f) 


a'(t)/(0). 


(30) 


The  prime  '  denotes  differentiation.  The  above  equations  are  sufficiently  general  for  computing  either 
the  off-time  or  on-time  response  for  an  arbitrary  transmitter  waveform.  Here  the  terms  off-  and  on- 
time  represent  measurement  times  taken  when  the  transmitter  current  source  is  either  on  or  off.  The 
pulse  I  (f )  function  we  consider  can  be  rectangular,  trapezoidal,  or  mixed  with  an  exponential  rise  time 
and  a  linear  ramp-off  time. 

To  begin,  we  consider  a  simple  case  where  the  field  response  can  be  described  by  an  analytical 
expression.  That  is,  a  VMD  embedded  in  a  uniform  conducting  medium  is  excited  by  a  unit  rectangular 
current  function  /(f) 


m 


0  for  t  <  0, 

<1  for  0  <  f  <  f  w , 
0  for  t  >  tw 


(31) 


where  tw  is  the  pulse  width.  In  the  spherical  coordinate  system,  the  magnetic  field  response  at  a 
distance  r  from  the  VMD  in  6  direction  is  given  by  (Wait,  1953) 

He  (r,  t )  =  sin  QRe  (t) ,  (32) 

4?rrJ 
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Figure  13:  Differential  Responses  derived  from  those  in  Fig.  12.  For  Rxl-Rx5  pair:  (a)  x-component; 
(b)  y-component;  (c)  z-component.  For  Rx7-Rx9  pair:  (d)  x-component;  (e)  y-component;  (f)  z- 
component.  D15XX-tot  and  D15XX-sca  denote  the  differential  responses  between  the  scattered 
and  the  total  x-components  in  the  three-layered  medium.  The  same  notation  rule  applies  to  the  y-  and 
z-components. 


where  mo  is  the  dipole  moment  and  9  is  the  polar  angle.  The  transient  characteristics  of  the  field  are 
described  by  the  function  Rg  (t ) : 


Rg(t) 


Cg  (t)  for  0  <  t  <  t w 

Cg{t)  -  Cg(t  ~  tw)  for  t  >  tw, 


(33) 


where  Cg{t)  is  a  step-on  response  function  (Wait,  1953).  Eq.  (33)  shows  that  the  field  responses  for 

a  dipole  excited  by  a  rectangular  current  can  be  expressed  as  a  supeiposition  of  the  two  step-function 

2 

responses.  The  function  Cg(t)  with  [3  =  is  given  by 

Cg{t)  =  erfc(/3)  +  2ir~*(p  +  2/33)  exp(-/32),  (34) 
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where  erfc(/3)  is  the  complement  of  the  error  function  of  argument  3.  Thus  we  can  use  Rg(t)  to 
examine  the  behavior  of  an  electromagnetic  wave  pulse  propagating  in  a  conducting  medium. 

For  an  exciting  rectangular  pulse  width  of  tw  =  10  ms  in  Eq.  (31),  Fig.  14  presents  Rg(t)  for 
different  host  media  conductivities  of  a  =  0, 0.2, 1,  2, 4  S/m.  When  r  =  20  m,  we  can  see  from  Fig. 
14  (a)  that  the  shape  of  the  on-time  responses  for  a  conductive  medium  is  not  rectangular  initially. 
Rather  it  surges  to  exceed  the  source  current  value  and  then  gradually  falls  back  and  settles  down. 
As  the  conductivity  of  the  host  increases,  the  on-time  fields  become  more  spread  out  and  the  peak 
amplitude  response  is  pushed  to  later  in  time.  Fig.  14  (b)  shows  these  same  plots  with  a  logarithmic  x- 
axis.  For  the  off-time  responses  shown  in  Fig.  14  (c),  we  see  that  the  effects  of  increasing  conductivity 
delay  the  responses  and  sign  reversals  to  a  later  time.  These  analytical  off-time  results  are  consistent 
with  those  of  the  numerical  simulation  (Song  et  ah,  2016).  When  considering  a  smaller  observation 
distance  r  =  2  m  and  contrasting  with  Figs.  14  (a)  and  (b),  we  observe  that  in  Figs.  14  (d)  and  (e) 
that  the  on-time  field  disperses  in  a  narrow  early  time  range  and  quickly  approaches  that  of  the  source 
waveform.  For  a  near- field  observation,  the  decaying  off-time  responses  in  Fig.  14  (f)  also  shift  toward 


Time  (ms) 


(a) 


(b) 


(c) 


(g) 


Figure  14:  Field  dispersion  function  Rg(t)  w.r.t.  conductivities  for  tw  =  10  ms.  r  =  20  m:  (a)  Finear 
scale,  (b)  Semilogx  scale,  (c)  Off-time  response.  r  =  2m:  (d)  Finear  scale,  (e)  Semilogx  scale,  (f) 
Off-time  response,  (g)  Field  dispersion  function  R$(t)  w.r.t.  distances  for  conductivity  of  4  S/m  and 

tw  =  1  ms. 
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early  times  as  compared  with  Fig.  14  (c).  This  shows  that  propagation  effects  in  a  conducting  medium 
displays  similar  characteristics  to  the  stand-off  distance  effects  studied  in  Sec.  3.2  using  a  step-off 
waveform. 

For  a  short  rectangular  pulse  of  1  ms,  Fig.  14  (g)  shows  the  dispersion  with  different  observation 
distances  from  0.5  m  to  20  m  for  the  medium  conductivity  of  4  S/m.  These  results  indicate  that, 
for  a  given  conductivity,  the  field  dispersion  phenomena  can  become  significant  for  larger  observation 
distances.  But  for  our  typical  sensing  range  within  1  or  2  m,  the  on-time  dispersion  virtually  disappears 
by  the  end  of  the  on-time  for  pulse  width  as  short  as  1  ms.  This  means  the  on-time  field  dispersion  is 
unlikely  to  be  of  concern  when  dealing  with  rectangular  pulse  source  functions. 

Let  us  now  examine  an  more  realistic  pulse  waveform  comprising  an  exponentially  rising  on-time 
and  linear-ramp  off  defined  as: 


m 


0  for  t  <  0, 

1  —  e~t/T  for  0  <  t  <  tw, 

. 

l-^(t-tw)  for  tw  <  t  <  tw  +  4, 
0  for  t  >  tw  +  tb 


(35) 


where  r  and  4  are  a  time  constant  and  slope  that  control  how  fast  a  current  is  switched  on  and  off. 
With  small  t  and  4  values,  7(f)  in  Eq.  (35)  approaches  an  ideal  rectangular  pulse.  We  simulate  four 
different  source  waveforms  given  tw  =  10  ms:  (a)  rectangular;  (b)  r  =  0.5  ms,  tb  =  50/rs;  (c)  t  =  1 
ms,  tb  =  50/rs;  (d)  r  =  1  ms,  4  =  1  ms,  shown  in  the  top  row  of  Fig.  14.  For  a  VMD  transmitter 
in  a  uniform  conducting  background  of  4  S/m  and  an  observation  at  r  =  20  m,  we  compute  the  field 
with  Eq.  (29)  by  convolving  the  waveform  function  of  Eq.  (35)  with  the  step-on  response  in  Eq.  (34). 
For  a  rectangular  pulse,  the  numerical  result  (bottom)  in  Fig.  15  (a)  is  the  same  as  the  analytical  one 
in  Fig.  14  (a):  a  dispersion  peak  appears  in  the  on-time  field  and  overlaps  with  the  step-on  response 
(gray  curve).  However  when  a  pulse  is  turned  on  exponentially  with  r  =  0.5  ms,  the  on-time  field 
(bottom)  in  Fig.  15  (b)  ramps  up  to  a  small  peak.  Increasing  r  =  1  ms,  the  dispersion  peak  is  almost 
completely  smoothed  out.  On  the  other  hand,  the  off-time  responses  visually  show  little  changes  from 
4  =  50  ps  to  ti,  =  1  ms.  These  results  show  that  for  realistic  pulse-waveforms  the  on-time  field 
dispersion  and  peak  are  greatly  diminished. 

Next  we  focus  on  the  off-time  responses  when  varying  the  pulse  width.  For  a  uniform  conducting 
medium  of  4  S/m,  we  compute  the  off-time  fields  Hg(t)  at  6  =  90°  and  at  different  observation  dis¬ 
tance  for  three  pulse  widths  of  tw  =  1  ms,  tw  =  10  ms,  and  tw  =  40  ms.  The  results  are  shown  in  Figs. 
16  (a)-(c).  For  a  given  pulse,  the  off-time  fields  sign  reversal  observed  at  each  location  occur  at  about 
the  same  time  for  each  pulse-width.  Also,  the  off-time  responses  at  different  observations  locations 
(0.5  m  -  20  m)  are  approximately  the  same  after  a  couple  of  milliseconds.  Again  the  analytical  field 
behavior  is  similar  to  the  numerical  study  of  offset  effects  on  measured  z-components  presented  in 
Sec. 3.2.  The  most  significant  effect  of  the  changing  pulse  width  is  the  change  is  the  change  in  the  field 
decay  rate.  The  decay  rate  t-5/2  of  an  impulse  response  is  shown  in  each  subplot,  and  it  is  evident 
that  the  decay  of  the  narrow  pulse,  tw  =  1  ms,  follows  this  same  relationship.  This  is  to  be  expected 
since  a  narrow  pulse  more  closely  resembles  an  impulse  source  excitation  that  do  the  longer  pulses. 
When  increasing  the  pulse  width  to  tw  =  10  ms  and  tw  =  40  ms,  we  see  that  the  field  decay  becomes 
slower.  By  viewing  a  rectangular  pulse  as  a  summation  of  two  opposite  polarities  of  step-on  functions, 
we  can  interpret  the  field  as  a  superposition  due  to  two  excitations  of  different  polarity.  For  a  longer 
pulse  duration,  the  “destructive”  superposition  occurs  later  in  time.  As  a  result,  the  field  would  behave 
more  like  the  one  due  to  a  step-off  excitation. 
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Figure  15:  Effects  of  pulse  shape  on  field  dispersion  at  r  =  20  m  in  a  uniform  medium  of  a  =  4  S/m 
with  excitation  waveforms  on  the  top  and  fields  on  the  bottom  row.  (a)  Rectangular  pulse,  (b)  r  =  0.5 
ms,  tf,  =  50  /ts;  (c)  r  =  1  ms,  tb  =  50  /is;  (d)  r  =  1  ms,  4  =  1  ms. 


(a)  (b)  (c) 


Figure  16:  Off-time  response  Hg(t)  w.r.t.  distances  and  pulse  widths  for  a  uniform  conducting 
medium  of  4  S/m.  (a)  tw  =  1  ms.  (b)  tw  =  10  ms.  (c)  tw  =  40  ms. 


Next  we  study  the  effects  of  pulse  width  on  the  background  and  scattered  fields  measured  in  a 
marine  environment  using  a  rectangular  Tx  loop.  As  before,  the  conductivities  of  the  seawater  and 
sediment  are  4  S/m  and  0.4  S/m.  The  sea  depth  is  6  m  and  sensor  height  is  1  m  above  the  sea  floor. 
For  a  z-directed  Tx  loop  of  2  m  x  1  m  and  a  centered  Rx  coil  of  0.15  x  0.15  m,  we  compute  the  fields 
with  the  two  rectangular  pulse  widths  of  tw  =  1  ms  and  tw  =  10  ms  by  the  convolution  of  Eq.  (30). 
Fig.  17  shows  the  background  and  scattered  fields  including  the  step-off  excitation.  Similar  to  the 
VMD  in  an  uniform  conducting  medium,  the  longer  duration  of  a  pulse  yields  strong  or  slow  decay 
responses.  It  also  shows  that  the  responses  with  tw  =  10  ms  better  approximate  those  of  an  ideal 
step-off  function  and  that  a  longer  pulse  enhances  the  target  response. 
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Figure  17:  Background  and  scattered  responses  with  a  Tx  loop  of  2  m  x  1  m  and  a  centered  Rx  coil  of 
0.15  m  x  0.15  m.  ttx  =  r rx  =  (0.0,  0.0,  dsea  —  h)  m,  h  =  1  m.  The  three-layer  parameters:  o\  =  0 
S/m,  02  =  4  S/m,  0-3  =  0.4  S/m,  dsea  =  6  111.  The  parameters  of  a  sphere:  r  =  (0.05,  0.05,  dsea  + 
0.30)  m,  a  =  0.05  m,  a  =  3.5  x  107  S/m,  \ir  =  300. 

3.6  Inversion  of  data  collected  underwater 

We  have  learnt  how  the  duration  of  a  pulse  can  affect  the  scattered  field  from  a  target  with  a  simulated 
2  111  x  1  111  Tx  loop  case.  In  this  section,  we  re-examine  this  phenomena  by  assuming  a  marine 
TEMTADS  2x2  sensor  system  (TT22)  excited  by  ideal  step-off,  rectangular  pulses  with  tw  =  1  ms 
and  tw  =  10  ms,  respectively.  Our  goal  is  to  evaluate  the  effectiveness  of  inversions  using  synthetic 
marine  scattered  data  when  a  magnetic  dipole  model  in  free  space  is  used  as  a  model  of  the  scattered 
field  response. 

A  TT22  system  shown  in  Fig.  18  consists  of  four  ^-directed  transmitters  of  dimension  0.35  111  x 
0.35  m  and  four  triaxial  receiver  cubes  of  0.1  111  x  0.1  111  x  0.1  m.  We  assumed  the  TT22  is  immersed 
in  sea  water  of  dsea  =  20  m  and  is  set  above  the  seafloor  at  a  height  of  0.2  m  to  interrogate  a  sphere 
buried  in  the  sediment.  The  sphere  has  a  radius  of  a  =  0.03  111  with  a  =  3.5  x  107  S/m,  /xr  =  300. 
Its  center  is  at  (0.05,  0.05,  dsea  +  0.30)  m.  The  conductivities  of  the  three  layers  are  a\  =  0  S/m  (air), 
0-2  =  4  S/m  (sea),  0-3  =  0.4  S/m  (sediment). 

With  the  different  driving  source  functions,  Fig.  19  (a)  shows  the  scattered  transient  responses 
measured  from  all  12  receivers  after  firing  one  transmitter,  e.g.,  Txl.  We  see  that  the  scattered  re¬ 
sponses  from  the  pulse  width  of  10  ms  are  stronger  than  those  from  the  pulse  width  of  1  ms,  particu¬ 
larly  for  late  times. 

Next  we  consider  a  simple  inversion  test.  For  terrestrial  UXO  sensing,  an  equivalent  induced 
dipole  model  has  been  used  effectively  to  describe  low  frequency  EMI  scattering  of  a  metal  target 
(Das  et  ah,  1990;  Pasion  and  Oldenburg,  2001;  Bell  et  ah,  2001;  Smith  et  ah,  2004).  In  the  time 
domain,  a  target  is  mathematically  characterized  by  a  3  x  3  magnetic  polarizability  tensor  P  and  has 
an  eigen-decomposition 


3 


(36) 


i= 1 
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where  &i{i  =  1,  2,  3)  is  the  orthonormal  eigenvector  representing  the  ith  principal  direction  of  dipolar 
polarization  with  respect  to  a  reference  system  and  Lj(i)  is  the  principal  polarization  strength  that  is 
a  function  of  the  geometry  and  material  properties  of  a  target.  The  task  of  EMI  sensing  of  UXO  is  to 
recover  the  principal  polarizabilities  from  measurements. 

To  solve  for  the  parameters  in  Eq.  (36)  the  following  system  of  equations  (Song  et  ah,  201 1)  are 
constructed 

d(i)  =  A(r,rs)q(f),  (37) 

where  d (t)  =  [di(t),  •  •  •  ,  d<M{t)]T  is  an  M  x  1  measured  data  vector  at  time  instant  t,  and  A(r,  rs)  is 
a  M  x  6  matrix  denoting  the  sensitivities  of  the  M  sensors  at  rs  to  the  object  located  at  r,  and  q(f )  a 
6x1  column  vector  whose  components  are  the  elements  of  the  polarizability  tensor  P(t)  of  an  object. 
They  are  given  by 
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where  aj(r,  rRxi,  rTXi)  is  6  x  1  column  vector  representing  spatial  sensitivities  of  the  zth  pair  of 
transmitter-receiver  to  the  object  located  at  r.  Notice  that  the  primary  fields  Hy  =  [HR,  Hl'r .  HR\  and 

Hr  =  [HR,  HnR.  H R]  in  Eq.  (38)  are  static  and  computed  for  a  transmitter  and  receiver  loop  in  free 
space. 

We  are  interested  to  know  if  Eqs.  (37)  and  (38)  be  directly  used  to  invert  the  scattered  data  in  the 
marine  environment  without  modification.  With  the  simulated  marine  TT22  in  Fig.  18,  we  proceed  to 
inversion  tests  by  using  the  three  sets  of  synthetic  data  in  Fig  19  (a):  the  step-off  and  two  rectangular 
pulses  with  widths  of  1  ms  and  10  ms.  At  this  stage  the  test  is  simplified  to  a  linear  inverse  problem  by 
assuming  a  known  location  of  the  target.  Fig  19  (b)  shows  the  recovered  polarizabilities  for  the  three 
excitation  waveforms.  For  each  current  function,  the  inversion  derives  three  identical  polarizabilties 
that  cannot  be  discerned  in  the  plot.  This  is  exactly  what  we  expected.  In  addition,  the  impact  of  pulse 
width  is  exhibited  in  the  recovered  polarizabilities.  The  longer  duration  of  an  excitation,  e.g.,  10  ms, 
produces  a  slower  decay  polarization  in  late  times  than  that  with  1  ms.  The  larger  width  of  pulse,  the 
better  the  recovered  polarizabilities  approach  those  from  the  step-off  driving  function. 

The  preliminary  inversion  test  illustrate  that  a  magnetic  dipole  model  in  free  space  is  readily 
adapted  to  processing  the  off-time  scattered  data  of  a  highly  conducting  and  permeable  object  buried 
in  the  sea-bottom  with  EMI  measurements  collected  in  a  conducting  underwater  environment. 
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Figure  18:  A  simulated  marine  TT22  system  consisting  of  four  planar  transmitter  (each  size  of  0.35 
m  x  0.35  m)  and  four  tri-axial  receiving  cubes  (each  size  of  0.10  m  x  0.10  m  x  0.10  m).  Sensor 
height  h  =  0.2  m.  The  center  of  the  sensor  is  rs  =  (0.0,  0.0,  dsea  —  h)  m.  The  three-layer  parameters: 
<7 1  =  0  S/m,  (72  =  4  S/m,  cr\  =  0.4  S/m,  dsea  =  20  m.  The  parameters  of  a  sphere:  a  =  0.03  m, 
a  =  3.5  x  10'  S/m,  nr  =  300.  Its  center  is  r  =  (0.05,  0.05,  dsea  +  0.30)  m. 
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Figure  19:  Given  the  setup  in  Fig.  18:  (a)  scattered  transient  responses  under  rectangular  pulse  width 
of  1  ms,  10  ms,  and  step-off  excitation  waveforms;  (b)  recovered  polarizabilities. 
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4  Numerical  investigation  of  an  insulated  VMD  immersed  in  a  uniform 
conducting  medium 


So  far,  our  studies  on  the  transient  background  and  scattered  responses  have  assumed  that  the  dipole  or 
loop  source  and  receiver  are  contained  in  a  uniform  or  layered  conducting  medium.  In  practice,  these 
sources  and  receivers  are  waterproofed  and  have  no  direct  contact  with  the  medium  under  study.  The 
fields  excited  by  an  insulated  loop  antenna  may  be  different  from  those  of  a  similar  bare  loop  antenna. 
In  this  section,  we  investigate  the  problem  by  enclosing  a  VMD  transmitter  in  a  spherical  cavity  that 
is  immersed  in  a  uniform  conducting  medium. 

Consider  a  magnetic  dipole  situated  at  the  origin  and  oriented  in  the  polar  direction  in  a  spherical 
co-ordinate  system  of  ( r,8,q i>)  in  Fig.  20.  According  to  Wait  (1952),  the  magnetic  dipole  in  this 
situation  will  generate  transverse-electric  spherical  waves  inside  and  outside  the  spherical  cavity.  For 
the  background  and  cavity  with  the  electrical  parameters  of  (ei,  01,  // 1 )  and  (e2, 02,  M 2),  respectively, 
the  resulting  non- zero  component  fields  at  r  in  the  frequency  domain  are 


Hlr  =  cos  0(1  +  7P")  exp(-7;r) 

Hie  =  ^3  sin  6(1  +  7 tr  +  7 fr2)  exp(~7 tr) 

TTli 

Ei,/,  =  sin 0(1  +  7;r)exp(-7;r), 


(39) 


where  l  =  1,2  and  7;  =  [707^(07  +  jwe/)]1/2.  For  l  =  1,  Eq.  (39)  describes  the  primary  fields  of 
the  magnetic  dipole  inside  the  spherical  cavity  with  mi  =  Id  A  with  a  current  of  I  and  the  differential 


Figure  20:  Schematic  setup  for  a  VMD  insulated  by  a  spherical  cavity  that  is  immersed  in  an  uniform 
conducting  medium.  The  radius  of  the  cavity  is  b.  The  radial  distances  of  77  <  b  and  77  >  b  at  a  polar 
angle  of  6  denote  the  observations  inside  and  outside  the  cavity. 
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area  dA.  For  l  =  2,  Eq.  (39)  describes  the  fields  in  the  region  exterior  to  the  cavity  with  an  equivalent 
dipole  moment  m2  contained  within  the  spherical  cavity.  The  ratio  between  m2  to  m\  is  given  by 


m2  _  71 B 
mi  72  C' 


The  secondary  fields  inside  the  cavity  due  to  the  insulation  boundary  are  given  by 

HirtS  =  ~r~~  cos6,[f,1'(7ir)  -  hirfir)]^ 

JUH 1 

H\e,s  =  — -  sin Oi'iilir) 

jw/xi  r 

A 

E\^s  =  sin  (71T), 


(40) 


(41) 


where  %\  is  a  spherical  Bessel  function  of  order  one  of  the  first  type,  defined  as  i\  (z)  =  cosh  z  — 
sinh  z/ z.  The  primed  quantities  indicate  a  differentiation  with  respect  to  the  argument  of  the  func¬ 
tion.  The  coefficients  A  and  B  in  above  equations  are  determined  by  imposing  the  continuities  of  the 
tangential  electric  and  magnetic  fields  E&  and  Hg  across  the  boundary  at  r  =  b: 


B  =  c  71  [ki  (71  fry  (71  ft)  -  4(71^1(7^)] 

71^1(726)^(716)  -  72 Ki (726)*1  (716)] 

4  _  -6^1(726)  -  CKiirgb)  (42) 

*1(716) 

C  =  jw/rimi7i/47r, 


where  k\  (z)  =  (1  +  z)  exp(— z)/z  is  the  spherical  Bessel  function  of  order  one  of  the  third  type. 

For  a  spherical  cavity  with  a  radius  of  b  =  20  cm  immersed  in  a  uniform  conducting  medium 
of  4  S/m,  Fig.  21  shows  the  VMD-generated  fields  at  r  =  b  and  6  =  45°  satisfying  the  boundary 
conditions:  E\^s  +  Ei^  =  and  H\g^s  +  H\g  =  ff2e»  in  frequency  domain,  which  involve  Eqs. 
(39),  (41),  and  (42). 

To  investigate  the  effect  of  the  cavity  we  compute  the  ratio  of  m2 /mi  using  Eq.  (40)  across  a 
wide  frequency  range  (Fig.  22).  For  a  20  cm  cavity  (Fig.  22a)  the  real  part  of  the  ratio  at  b  =  20 
cm  has  the  constant  value  of  1  across  a  wide  range  of  frequencies  and  the  imaginary  part  is  less  than 
1,  linearly  increasing  with  logarithmic  frequencies.  The  ratio  starts  to  oscillate  at  high  frequencies 
>  10  MHz.  When  decreasing  the  size  of  the  cavity  to  b  =  10  cm  and  6  =  5  cm,  the  oscillation 
in  the  ratio  moves  toward  higher  frequencies  after  100  MHz  and  the  amplitude  becomes  smaller.  At 
sufficiently  low  frequencies,  the  equivalent  dipole  moment  m2  for  an  insulated  VMD  approaches  mi 
for  an  uninsulated  VMD.  We  conclude  that  the  cavity  has  little  effect  on  observed  fields  in  the  exterior 
region. 

We  now  turn  our  attention  to  the  fields  inside  the  cavity.  Referring  to  Fig.  20,  for  6  =  0.2  m  and 
0  =  45°  we  compute  the  transient  fields  inside  the  cavity  at  r\  =0.1  m  and  the  fields  outside  the  cavity 
at  r2  =  1  m,  r2  =  5  m  and  r2  =  10  m.  Fig.  23  presents  vertical  magnetic  fields  inside  and  outside 
of  the  cavity.  The  fields  outside  with  and  without  the  cavity  are  the  same,  consistent  with  our  analysis 
in  the  frequency  domain.  For  the  fields  inside,  we  plot  the  primary  and  scattered  field  separately  for 
better  visualization  of  their  behavior.  The  results  shows  that  the  primary  field  inside  is  far  weaker  than 
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Figure  2 1 :  A  check  of  boundary  conditions  of  an  insulated  VMD  in  frequency-domain  for  9  =  45° 
and  b  =  0.2  m.  (a )E(/).  (b )Hq. 


the  scattered  field  inside.  Surprisingly,  the  scattered  field  inside  at  r\  =0.1  m  (thick  curve)  almost 
overlaps  with  the  field  outside  at  =  1  m  roughly  beginning  with  0.01  ms.  Increasing  r2  to  5  m  and 
10  m,  we  see  that  it  takes  longer  for  the  fields  inside  and  outside  to  overlap  with  convergence  occurring 
somewhere  between  0.5  ms  and  1  ms.  This  correspondence  of  background  fields  across  widely  offset 
observation  distances  is  similar  to  what  we  have  observed  in  previous  sections  when  numerically  and 
analytically  studying  the  offset  effects  of  bare  loops  or  VMD  in  layered  or  uniform  conducting  media. 
The  cavity  insulation  appears  to  have  little  effect  on  the  vertical  fields  both  inside  and  outside  the 
cavity. 
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Figure  22:  Ratio  of  dipole  moments  m2 /mi  for  an  insulated  VMD.  The  radius  b  of  a  spherical  cavity: 
(a)  20  cm;  (b)  10  cm;  (c)  5  cm. 
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(a)  (b)  (c) 

Figure  23:  cLHz/dt  responses  inside  and  outside  a  spherical  cavity  defined  in  Fig.  20.  With  b  =  0.2 
m,  0  =  45°,  and  n  =  0.1  m,  r2  =  (a)  1  m;  (b)  5  m.  (c)  10  m.  Field  inside  =  Primary  field  inside  + 
Scattered  field  inside. 
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5  Underwater  measurements 


The  original  intent  of  this  project  was  to  utilize  data  collected  by  marine  sensors  under  development 
as  part  of  other  funded  SERDP/ESTCP  projects.  However,  those  projects  have  been  delayed  and 
suitable  data  was  not  available.  Through  Black  Tusks  interactions  with  Gap  Explosive  Ordnance 
Detection  (GapEOD)  in  Brisbane,  Australia,  an  opportunity  arose  to  collect  some  marine  data  with  an 
EM  system. 

For  the  marine  measurements,  GapEOD  built  a2mxlmxlm  fiberglass  frame  and  encased  two 
receiver  cubes  in  epoxy  to  make  them  waterproof  (Fig.  24.)  The  transmitter  loop  comprised  12  turns 
of  wire  arranged  in  a  2  m  by  1  m  rectangle  at  approximately  the  same  vertical  level  as  the  base  of  the 
receiver  cubes.  A  24V  power-supply  was  used  to  provide  a  maximum  current  of  1 1.4  A  (at  23.4  V)  to 
a  Gap  LPTX-25  transmitter  that  operated  at  a  25  Hz  base-frequency  with  a  50%  duty  cycle  waveform. 
Four  separate  measurements  campaigns  were  conducted  either  within  Moreton  Bay,  close  to  Brisbane 
or  the  Gold-Coast  Seaway  near  Surfers  Paradise: 

1)  Measurements  over  a  range  of  ordnance  items  between  water  depths  of  2.5  and  14  m  conducted 
over  two  days  (13  and  16  March  2015); 

2)  Terrestrial  and  marine  measurements  conducted  over  a  105  mm  projectile  at  exactly  the  same 
locations  relative  to  the  sensor  frame  (July  2015). 

3)  More  controlled  measurements  over  a  105  mm  projectile  in  shallow  water  (April  2016); 

4)  Terrestrial  and  marine  coil  and  B-field  full-waveform  measurements  to  better  understand  the 
underwater  noise  environment  and  primary  field  distortions  caused  by  the  conductive  sea-water 
(June  2016). 

5.1  Measurement  campaign  1:  Background  and  proof  of  concept 

For  the  first  measurement  campaign  we  chartered  a  small  vessel  (Fig.  25)  and  collected  background 
data  at  various  different  water  depths  (Fig.  26)  within  Moreton  Bay.  Water  samples  taken  during  the 


Figure  24:  Frame  built  for  the  underwater  tests,  comprising  two  receive  cubes  and  a  2  m  x  1  m 
transmitter  coil. 
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UltraTEM  measurements  revealed  that  the  sea-water  had  a  conductivity  of  6  S/m.  For  the  measure¬ 
ments  the  two  receivers  were  placed  80  cm  apart  and  equally  offset  by  40  cm  from  the  transmitter  loop 
center.  The  background  signal  in  the  horizontal  components  is  quite  weak:  this  is  because  at  times 
greater  than  100  ps,  the  background  behaves  as  if  it  were  generated  by  a  propagating  smoke-ring 
(Nabighian,  1979)  whose  dimensions  are  several  times  larger  than  the  size  of  the  transmitter.  The 
field  near  the  center  of  that  smoke  ring  is  predominantly  vertical.  The  observations  show  that  the 
vertical  component  field  decays  relatively  rapidly  with  time:  at  a  rate  of  between  t-5/2  and  t-3.  The 
amplitude  of  the  decay  increases  with  water-depth  but  converges  to  a  steady  state  value  equivalent 
to  an  infinite  half-space.  These  decay  rates  agree  with  theoretical  predictions  and  modeling  results 
developed  within  this  project  that  account  for  sea-water  interface  effects  (section  3.1).  The  interface 
effects  account  for  the  slightly  faster  decay  rate  of  the  measurements  collected  in  shallower  water. 

We  also  collected  data  with  ordnance  items  placed  inside  the  PVC  tube  with  the  sensor  frame 
positioned  on  the  sea-bottom.  For  those  measurements  we  first  collected  a  background  measurement, 
then  raised  the  frame  to  the  surface  and  inserted  the  ordnance  item  or  surrogate  into  the  PVC  tube, 
before  the  frame  was  lowered  over  approximately  the  same  location.  Measurements  were  then  taken 
with  the  frame  stationary  for  at  least  a  minute  per  sample  location.  Photos  of  some  of  the  items 
measured  during  Measurement  Campaign  1  are  shown  in  Fig.  27  with  the  resulting  receiver  data 
shown  in  Fig.  28  .  All  these  measurements  were  collected  when  the  ordnance  item  or  surrogate  was 
at  either  end  of  the  PVC  tube  (Pos.  1  directly  under  the  transmitter  wire,  or  Pos  2.  1  m  outside  the 
transmitter  frame).  In  each  case,  we  plot  the  data  from  the  closest  receiver  which  is  still  always  at 
least  1  m  away  from  the  object  being  measured.  At  early  times  the  metal  and  background  signals  are 
comparable,  but  the  sea-water  response  dies  down  much  faster  than  the  response  from  the  metallic 
items.  Hence,  good  strong  signals  from  the  metallic  items  are  evident  at  later  times.  These  results 
demonstrate  that  detection  of  metallic  ordnance  in  sea-water  is  possible  using  a  large  transmitter  loop 
and  three-component  receiver  cubes. 
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Figure  25:  Underwater  tests  were  conducted  in  Moreton  Bay  near  Brisbane,  (a):  System  being  moved 
over  the  side  of  the  vessel  prior  to  lowering  into  the  water;  (b)  Components  of  the  EM  system  including 
power-supply,  LPTX-25  transmitter,  Toughpad  for  display  of  data  and  UltraTEM  receiver  system. 
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Figure  26:  Background  signals  at  different  water  depths  as  measured  by  the  UltraTEM  system  (left 
X-component,  center  Y-component  and  right  Z-component). 


Figure  27:  Items  measured  during  phase  1,  including  81  mm  mortar  (top  left),  ferrous  cylinder  (top 
right)  and  105  mm  projectile  (bottom  left).  An  image  of  the  sample  tube  with  a  105  mm  projectile 
present  is  at  the  bottom  right. 
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(b)  105  mm  projectile 


Figure  28:  Underwater  measurements  over  several  test  items  at  Pos  1  (under  the  edge  of  transmitter 
coil)  for  (a)  and  first  measurement  in  (b)  with  second  measurement  at  Pos  2(1  m  from  edge  of  the 
transmitter  coil  and  outside  the  loop). 


SERDP  2412  Final  Report 


40 


Nov-2016 


5.2  Measurement  campaign  2:  Comparison  of  terrestrial  and  marine  measurements 


During  measurement  campaign  2  we  collected  data  with  exactly  the  same  transmitter/receiver/object 
geometry  in  both  the  terrestrial  and  marine  environments  (Fig.  29).  All  measurements  were  collected 
over  a  105  mm  projectile.  At  later  times  (beyond  0.5  to  1  ms)  the  terrestrial  and  marine  results  are 
virtually  identical  (Figs.  30  &  31).  This  means  that  from  a  modeling  perspective,  these  later  time 
data  can  be  treated  the  same  in  both  the  terrestrial  and  marine  environments,  which  agrees  with  the 
conclusions  from  our  modelling  study  and  previous  SERDP  sponsored  research.  The  measurements 
with  the  105  mm  in  the  marine  environment  undershoot  the  background  measurements  at  early  time 
(Fig.  31(a)),  meaning  that  background  subtraction  reveals  negative  transients  from  the  projectiles 
(Fig.  31(a)).  This  interaction  phenomenon  has  been  observed  previously  (e.g.  Shubitidze,  2011)  and 
increases  with  increasing  offset  between  the  object  and  receiver.  Flowever,  the  interaction  effects 
evident  in  these  measurements  are  much  later  in  time  than  have  been  observed  previously.  We  suspect 
these  are  due  to  incorrect  background  measurement  caused,  in  part,  by  the  requirement  to  lift  and  then 
reposition  the  frame  between  measurements. 


Figure  29:  Layout  of  the  frame  used  for  the  second  measurement  campaign  showing  the  locations  of 
measurements  collected  over  a  105  mm  projectile  in  both  the  terrestrial  and  marine  environments. 


SERDP  2412  Final  Report 


41 


Nov-2016 


X 


Y 


Z 


Figure  30:  Terrestrial  measurements  of  a  105  mm  projectile  at 
sensing  frame. 


three  different  positions  under  the 
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Figure  31:  Marine  and  terrestrial  measurements  are  virtually  identical  from  about  0.5  to  1  ms  and 
later  in  time,  with  the  seawater  effects  predominantly  occurring  at  early  time.  We  suspect  that  the 
undershoots  that  occur  at  early  time  are  due  to  poor  background  removal  and  are  not  indicative  of  the 
true  response  from  the  object. 
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5.3  Measurement  campaign  3:  Early  time  effects 


In  the  third  measurement  campaign  we  focused  on  understanding  the  early-time  effects  that  were 
apparent  in  the  second  measurement  campaign  but  which  we  suspected  were  due  to  poor  background 
removal.  The  transmitter  and  receiver  set-up  was  the  same  as  the  previous  two  campaigns  but  this 
time  we  installed  the  receiver  cubes  and  transmitter  loop  onto  a  smaller,  more  portable  fiberglass 
frame  (Fig.  32).  In  addition,  rather  than  retrieving  the  frame  between  measurements,  we  kept  the 
frame  in  a  fixed  location  and  used  a  diver  to  move  the  ordnance  item  (a  105  mm  projectile)  to  the 
different  measurements  positions  shown  in  Fig.  33.  Between  each  measurement  on  the  105  mm 
projectile  we  collected  a  background  measurement  so  were  better  able  to  remove  backgrounds  than 
we  were  in  the  second  campaign.  One  receiver  cube  was  placed  in  the  center  of  the  transmitter,  and 
the  second  was  placed  outside  the  transmitter  footprint,  2.0  m  from  the  center  of  the  system  on  the  left 
side  of  the  system.  The  ordnance  item  was  moved  between  three  different  lateral  position,  0.505  m 
right  of  the  center,  0.495  m  left  of  the  center  and  then  1.5  m  left  of  the  center  (0.5  m  from  the  second 
receiver  cube).  Measurements  were  conducted  in  3.3  to  4  m  water-depth  (the  tidal  variation  during  the 
measurement  period)  with  the  frame  on  the  sea-bottom  and  then  with  the  frame  raised  50  cm  above 
the  sea-bottom,  but  the  ordnance  item  kept  on  the  sea-bottom.  Water  salinity  was  again  quite  high 
with  a  measured  conductivity  reading  of  6.16  S/m  (similar  to  the  values  observed  in  Moreton  Bay  in 
the  previous  measurement  campaigns). 

The  response  from  cube  1  for  the  105  mm  projectile  in  position  2  when  the  system  is  elevated  50 
cm  above  the  sea- floor,  shows  almost  identical  responses  in  X-  and  Z-component  data  for  the  105  mm 
projectile  insulated  (inside  a  thick  plastic  bag)  and  not-insulated3  (Fig.  34).  The  Y-component  data 
are  different  but  we  suspect  that  is  due  to  slight  orientation  differences  between  the  insulated  and  non- 
insulated  samples.  The  background  corrected  responses  are  indistinguishable  from  what  they  would 
be  for  the  equivalent  terrestrial  measurement.  Any  early  time  effects  of  the  sea-water  (like  CCR) 
would  be  more  pronounced  for  larger  receiver  to  object  offsets.  Fig.  35  and  36  show  the  responses 
of  cubes  1  position  1  and  cube  2  position  3  for  the  insulated  and  non-insulated  105  mm  projectiles. 
These  measurement  scenarios  correspond  to  increased  receiver  to  object  offset,  in  one  case  with  the 

’The  projectile  has  a  layer  of  paint  so  is  already  partly  insulated 


Figure  32:  Test  frame  used  for  the  third  measurement  campaign.  A  diver  placed  the  105  mm  projec¬ 
tile  at  each  of  the  test  locations  and  background  measurements  were  collected  before  and  after  each 
measurement. 
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Figure  33:  Plan  view  of  the  system  and  source  locations  for  the  third  measurement  campaign. 

object  outside  the  transmitter  and  receiver  inside  and  the  other  with  the  reverse.  Insulated  and  non- 
insulated  responses  are  again  virtually  identical  to  each  other  and  to  the  equivalent  terrestrial  response 
with  small  deviations  at  early  time.  These  deviations  are  likely  due  to  imperfect  background  removal 
(the  background  is  more  than  2  decades  larger  than  the  UXO  response  at  early  time). 

Based  on  the  measurements  obtained  in  measurement  campaign  3  we  conclude  that  background 
corrected  terrestrial  and  marine  measurements  are  indistinguishable  at  times  less  than  200  microsec¬ 
onds  (and  probably  even  earlier),  and  the  105  mm  projectile  with  paint  covering  does  not  exhibit  any 
significant  CCR. 

As  part  of  the  third  measurement  campaign  we  also  collected  several  minutes  of  data  with  the 
transmitter  on  and  the  frame  on  the  beach  and  again  when  it  was  underwater.  Fig.  37  compares  the 
sensor  noise  of  the  underwater  and  terrestrial  measurements.  The  noise  levels  are  comparable  with 
a  very  slight  increase  in  early  time  noise  for  the  Z-component  receiver  data  when  the  system  was 
underwater. 
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Figure  34:  Comparison  of  insulated  and  non-insulated  response  for  position  2,  cube  number  1  (cube 
and  object  close  together)  with  the  system  elevated  above  sea-floor.  Top  plots  are  raw  data  and  bottom 
plots  are  after  background  correction.  The  mismatch  in  Y-components  is  probably  due  to  orientation 
differences  between  the  two  sets  of  measurements. 
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Figure  35:  Comparison  of  insulated  and  non-insulated  response  for  position  1,  cube  number  1  (cube 
and  object  medium  distance  apart)  with  the  system  on  the  sea-bottom.  Top  plots  are  raw  data  and 
bottom  plots  are  after  background  correction. 
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Figure  36:  Comparison  of  insulated  and  non-insulated  response  for  position  3,  cube  number  2  (cube 
and  object  far  apartjwith  the  system  on  the  sea-bottom.  Top  plots  are  raw  data  and  bottom  plots  are 
after  background  correction. 


Figure  37:  Comparison  of  sensor  noise  with  the  system  in  air  on  the  beach  and  underwater:  from  left 
to  right  are  the  X,  Y  and  then  Z-component  data. 
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5.4  Measurement  campaign  4:  Noise  and  primary  field  distortion 


For  the  forth  measurement  campaign  we  returned  to  the  same  beach-side  environment  of  the  third 
campaign.  Our  objectives  with  these  measurements  were  to  (1)  better  understand  the  noise  spectra  in 
the  underwater  environment;  and  (2)  verify  the  primary  field  dispersion  calculations  derived  in  Sec. 
3.5.  We  used  the  same  transmitter  set-up  as  previous  campaigns  but  augmented  our  coil  measurements 
with  flux-gate  magnetometer  measurements.  Both  the  coil  and  flux-gate  data  were  recorded  with  an 
Electromagnetic  Imaging  Technology  (EMIT)  SMARTem24  receiver:  for  the  coil  the  sample  rate  was 
120  kHz  while  for  the  flux-gate  the  sample  rate  was  24  kHz.  We  used  a  Bartington  MAG-03MSSL70 
flux-gate  for  the  measurements  and  a  nulling  control  box  (to  null  the  Earth’s  magnetic  field)  built 
by  EMIT.  A  series  of  measurements  were  collected  on  land  and  in  water  at  increasing  lateral  offset 
distances  from  the  transmitter  loop. 

5.4.1  Noise  and  primary  field  oscillations 

One  motivation  for  collecting  the  noise  spectra  were  some  results  reported  in  Bell  et  al.  (2016)  where 
a  significant  increase  in  noise  was  reported  when  a  thinly  insulated  transmitter  loop  was  immersed  in 
salt-water.  In  simple  terms  the  transmitter  can  be  modelled  as  a  parallel  RLC  circuit  (Fig.  38)  with 
damping  factor 


(43) 


with  L  the  inductance,  C  the  capacitance  and  R  the  resistance.  When  (  >  1  the  response  of  the 
transmitter  is  overdamped  and  the  field  decays  exponentially  on  transmitter  shut-off.  The  resistor 
value  for  the  transmitter  is  typically  chosen  so  that  the  system  is  either  critically  damped  ((  =  1)  or 
overdamped  (£  >  1).  The  wire  used  for  our  transmitter  loop  comprised  12  turns  of  1  mm2  wire  with 
a  measured  inductance  of  630  pH  and  in-air  capacitance  of  200  pF.  For  critical  damping,  a  damping 
resistance  of  approximately  890  C  is  required. 


I 


Figure  38:  Parallel  RLC  Circuit. 


The  sea-water  increases  the  capacitance  of  the  loop  through  increased  turn-to-core  capacitance 
(stray  capacitance  between  the  transmitter  wires  and  the  core  material).  When  a  poorly  insulated 
transmitter  wire  that  is  critically  damped  is  immersed  in  sea-water  the  increase  in  C  causes  £  to 


decrease  to  less  than  one  and  the  system  becomes  under-damped.  When  the  transmitter  turns  off  it 
oscillates  at  the  natural  resonance  frequency  of  an  RLC  circuit  of 


(44) 
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Based  on  the  inductance  and  capacitance  values  measured  in  air,  the  resonant  frequency  of  our  trans¬ 
mitter  loop  is  approximately  450  kHz.  These  oscillations  can  be  avoided  either  by  insulating  the 
transmitter  loop  (as  found  by  Bell  et  ah,  2016)  or  by  increasing  the  damping  resistor  in  the  system.  To 
test  this  latter  assertion,  we  collected  measurements  with  a  coil  on  land  and  in-water  using  our  thinly 
insulated  transmitter  loop  with  and  without  a  damping  resistor.  4  For  each  case  studied  we  collected  a 
minute  of  data  with  the  receiver  coil  in  the  center  of  the  2  nr  x  1  m  transmitter  loop.  Fig.  39  shows  the 
first  200  ps  after  transmitter  turn-off  (at  10  ms).  For  both  the  on  land  and  in-water  measurements  with 
no-damping  resistor,  oscillations  are  evident  and  of  increased  magnitude  in  sea-water  (£  decreases 
with  increasing  C).  Unfortunately,  the  band-width  of  the  system  was  insufficient  to  accurately  resolve 
the  resonance  frequency  (the  oscillations  are  under-sampled  but  clearly  present).  When  the  damping 
resistor  is  reintroduced  into  the  system,  the  oscillations  disappear  for  both  the  on-land  and  in-water 
measurements  and  all  that  remains  is  the  self-response  from  the  transmitter  and/or  receiver  coil.  On- 
land  and  in-water  responses  are  almost  identical. 

With  the  damping  resistor  in  the  system,  we  collected  additional  coil  measurements  at  increasing 
lateral  offsets  from  the  transmitter  coil  center  (Fig.  40).  The  horizontal  component  response  increases 
as  the  coil  is  moved  towards  the  edge  of  the  transmitter  wire,  then  flips  sign  (as  does  the  vertical 
component  response)  once  measurements  are  collected  outside  the  coil.  The  off-time  responses  (at 
least  on  the  linear  scale  and  limited  time -range  shown)  for  the  on-land  and  in-water  measurements  are 

4The  case  of  no  damping  resistor  is  equivalent  to  R  — >  oo  and  £  — s ►  0  and  hence  underdamped. 


(a)  Horizontal  field 


Time  (ms) 

(b)  Vertical  field 

Figure  39:  Off-time  coil  response  in  the  horizontal  (a)  and  vertical  (b)  directions  for  a  transmitter  loop 
with  and  without  a  damping  resistor  on-land  and  in  sea-water. 


SERDP  2412  Final  Report 


50 


Nov-2016 


almost  identical  for  all  offset  locations. 


- Land  center  coil 

- Land  0.5  m 

Land  1 .5  m 

- Land  2.0  m 

—  Land  4.0  m 

- Marine  center  coil  . 
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(a)  Horizontal  field 


Figure  40:  Off-time  coil  response  at  different  horizontal  offsets  on  land  and  in  sea-water  in  the  hori¬ 
zontal  (a)  and  vertical  (b)  directions  for  a  transmitter  loop  with  damping  resistors. 

We  also  computed  estimates  of  the  power-spectral  density  (PSD)  for  the  loop  with  and  without 
a  damping  resistor  present,  for  the  system  on  land  and  in-water,  as  well  as  without  the  transmitter 
present  (Fig.  41).  With  the  transmitter  off,  the  marine  environment  is  lower  noise  than  the  adjacent 
terrestrial  environment.  The  50  Flz  power-line  noise  and  its  harmonics  are  clearly  evident  but  so  also 
are  a  number  of  other  lower  frequency  noise  sources  (of  around  15  Hz  and  harmonics).  For  the  vertical 
field,  the  PSD  as  a  function  of  frequency  is  almost  identical  for  on-land,  in-water,  damping  resistor 
present  or  absent.  However,  for  the  horizontal  component  there  is  a  marked  increase  in  spectral  energy 
in  both  the  on-land  and  in-water  measurements  with  no-damping  resistor,  consistent  with  the  results 
reported  in  Bell  et  al.  (2016). 

Finally,  we  computed  the  PSD  for  the  same  lateral  offset  positions  of  Fig.  40.  For  both  the 
horizontal  and  vertical-field  the  on-land  and  in-water  PSD  responses  at  each  offset  are  almost  identical. 
From  these  results  we  can  conclude  that  there  are  no  significant  differences  between  on-land  and 
in-water  sensing  from  the  perspective  of  the  frequency  content  of  the  noise  and  signal  (for  the  low- 
frequency  range  considered  here). 

5.4.2  Primary  field  dispersion 

In  Sec.  3.5  we  considered  the  dispersion  of  the  primary  field  waveform  in  the  conductive  seawater 
environment.  Our  modelling  of  an  ideal  step-off  waveform  in  Fig.  14  showed  a  significant  dispersion 
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(b)  Vertical  field 


Figure  4 1 :  Power  spectral  density  in  the  horizontal  (a)  and  vertical  (b)  directions  for  receiver  coil  at 
the  center  of  a  transmitter  loop  with  and  without  a  damping  resistor  on-land  and  in  sea-water.  Also 
shown  are  the  power  spectral  densities  for  the  coil  in  land  and  in  sea-water  with  no  transmitter  present. 

of  the  primary  field  with  increased  lateral  distance  from  the  transmitter.  This  effect  was  considerably 
diminished  when  accounting  for  the  exponential  rise  time  of  a  transmitter  loop  (Fig.  14(c)  &  (d)). 
To  confirm  these  theoretical  predictions  we  collected  one  minute  of  flux-gate  data  in  the  center  of  the 
loop  and  at  1  m  lateral  offset  distances  out  to  6  m  away  from  the  loop  both  underwater  and  on  land.  In 
all  cases  the  transmitter  was  operated  at  25  Hz  base-frequency  at  50%  duty  cycle  and  approximately 
2.6  A  of  current5.  The  flux-gate  dynamic  range  was  exceeded  for  the  center  loop  and  1  m  lateral  offset 
distances  so  these  were  excluded  from  the  plots  in  Fig.  43.  As  was  found  in  Sec.  3.5  the  primary 
field  dispersion  increases  with  increasing  lateral  offset  in  the  conductive  underwater  environment  but 
remains  virtually  unchanged  in  the  terrestrial  environment.  For  the  in-water  measurements,  even  at 
the  maximum  lateral  offset  distance  of  6  m,  the  effect  on  the  transmitter  on-time  is  marginal,  with 
the  largest  impact  occurring  during  the  rapid  shut-off  of  the  primary  field  at  the  end  of  the  transmit 
cycle.  In  most  UXO  detection  scenarios,  the  transmitter,  receiver  and  object  would  typically  be  within 
several  meters  of  each  other  and  we  conclude  that  the  primary  field  dispersion  is  not  an  important 
practical  consideration. 


5 


We  used  this  low  value  of  current  to  minimize  the  range  where  the  flux-gate  dynamic  range  was  exceeded 
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Figure  42:  Power  spectral  density  for  the  same  offset  locations  as  Fig  40. 
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Figure  43:  Normalized  magnetic  field  recorded  over  (a)  one  full  period  of  a  25  Hz  50%  duty  cycle 
transmitter  waveform  on  land  and  in  water  at  lateral  offset  distances  of  2,  3,  4,  5  and  6  m  both  under¬ 
water  (colored  lines)  and  on  land  (black  lines).  All  terrestrial  measurements  are  shown  with  the  same 
color  as  they  are  almost  identical,  (b)  is  a  close-up  of  the  on-time  waveform  and  (c)  shows  the  early 
off-time. 
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6  Conclusions 


In  this  report,  we  adapted  the  integral  equation  approach  to  calculate  EMI  responses  of  the  background 
and  scattered  objects  in  a  marine  setting.  In  the  formulation,  underwater  environments  are  treated  as  a 
three-layer  conducting  medium  that  consists  of  the  air,  the  sea,  and  the  ground.  A  highly  conducting 
and  permeable  sphere,  which  has  EMI  scattering  properties  expressed  analytically  in  the  frequency 
domain,  is  used  as  a  buried  target  in  sea  sediments.  The  transient  fields  are  obtained  by  digital  filtering 
techniques  over  the  harmonic  fields  for  either  dipole  or  loop  transmitters  and  receivers.  A  range  of  syn¬ 
thetic  experiments  were  conducted  with  respect  to  various  factors  that  might  influence  the  background 
and  scattered  responses,  including  the  sea  depth,  the  size  of  a  loop,  offset,  host  conductivity,  excita¬ 
tion  waveforms  and  receiver  and  transmitter  insulation.  These  simulation  studies  were  augmented  by 
a  number  of  measurements  using  a  custom  built  fibreglass  frame  with  a  2  m  x  1  m  transmitter  loop 
and  two  receiver  cubes. 

We  can  summarize  the  main  conclusions  of  this  report  as  follows: 

1)  Effect  of  sea-depth.  Simulation  studies  and  measurements  showed  that  the  background  response 
in  a  shallow  sea  decay  faster  than  in  a  deep  sea  (water  depth  50  m  or  more).  In  a  deep  sea 
where  the  sensor  is  far  away  from  the  air-sea  interface,  the  corresponding  background  responses 
asymptotically  approaches  the  response  of  a  half-space.  Observed  decay  rates  were  between 
t-5/2  and  t~3  depending  on  water  depth. 

2)  Effect  of  conductive  environment  on  scattered  response.  We  found  that  dispersion  effects  due 
to  conductivity  and  any  current  channel  response,  were  more  important  at  early  times  and  large 
receiver  to  object  distances.  However,  within  the  interesting  time  range  of  0.1  ms  -  25  ms,  the 
CCR  contributes  little  to  the  target  signals  and  thus  can  be  ignored.  At  times  beyond  several 
hundred  microseconds,  the  ECR  response  approaches  the  value  for  the  same  object  embedded 
in  free  space. 

3)  Transmitter  geometry  and  receiver  position.  The  background  responses  measured  inside  and 
outside  of  a  transmitter  loop  are  very  similar  and  at  later  times  vary  little  with  offset  distance 
from  the  transmitter  center.  A  large  transmitter  loop  can  excite  a  strong  scattered  response  from 
a  buried  object  but  also  induces  a  larger  background  response.  Due  the  small  variation  of  the 
background  response  with  position,  a  number  of  different  background  suppression  methods  can 
be  applied  to  make  measurements  with  large  loops  feasible.  Receivers  distributed  with  radial 
symmetry  around  the  transmitter  loop  experience  exactly  the  same  background  field  from  the 
conductive  sea-water.  Thus  differential  measurements  are  one  way  to  effectively  suppress  the 
response  of  the  seawater. 

4)  Field  dispersion.  Increasing  host  conductivity  and  offset  distances  increases  dispersion  of  the 
primary  field.  But  for  our  typical  sensing  range  within  1  or  2  m,  the  on-time  dispersion  virtually 
disappears  by  the  end  of  the  on-time  for  pulse  width  as  short  as  1  ms.  This  means  the  on- 
time  field  dispersion  is  unlikely  to  be  of  concern  for  practical  UXO  sensing.  Results  with  a 
simulated  marine  TEMTADS  system  showed  that  the  scattered  responses  excited  from  a  long 
pulse  of  10  ms  are  stronger  than  those  from  a  short  pulse  of  1  ms  in  particular  for  late  times. 
This  demonstrates  one  method  to  further  enhance  target  signals,  which  is  equally  effective  and 
equivalent  in  the  terrestrial  realm. 

5)  Inversion  tests.  Tests  with  simulated  marine  TEMTADS  data  showed  that  a  terrestrial  dipole 
polarizability  model  could  be  used  to  successfully  characterize  the  scattered  data  from  a  buried 
metallic  object  in  the  conducting  underwater  environment. 
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6)  Insulation  of  transmitter  and  receiver.  Our  preliminary  numerical  investigation  showed  that 
any  insulation  of  the  transmitter  or  receiver  will  have  negligible  effect  on  marine  UXO  sensing 
(assuming  the  transmitter  loop  has  been  adequately  damped). 

7)  Noise  and  spectral  measurements.  Underwater  measurements  with  no  transmitter  present  were 
lower  noise  across  most  of  the  frequency  spectrum  compared  to  the  on-land  measurements. 
When  the  transmitter  was  operating  and  sufficiently  damped  by  a  resistor,  the  spectrum  and 
off-time  behaviour  on  land  and  in  water  were  virtually  identical.  From  these  results  we  can 
conclude  that  there  are  no  significant  differences  between  on-land  and  in-water  sensing  from 
the  perspective  of  the  frequency  content  of  the  noise  and  signal. 

6.1  Future  research  directions 

Based  on  the  results  of  this  limited  scope  project,  we  conclude  that  underwater  detection  and  character¬ 
ization  of  buried  metal  using  EMI  based  sensor  is  feasible.  While  the  conductive  sea-water  introduces 
complications  to  the  measurement  process,  there  are  numerous  practical  options  for  mitigation  of  these 
effects.  The  conductive  sea-water  can  impact  the  scattered  field  from  a  buried  metallic  object,  but  typ¬ 
ically  only  at  very  early  times  and  at  large  receiver  to  object  offsets.  Neither  of  these  conditions  are 
commonly  encountered  in  practice.  Thus  apart  from  the  practical  considerations  of  operating  under¬ 
water,  the  modelling  techniques  and  methods  that  have  been  successfully  demonstrated  in  terrestrial 
environments  can  be  utilized  for  marine  detection  and  characterization. 

Once  we  had  obtained  an  understanding  of  the  impact  of  the  marine  environment  on  detection  and 
characterization  of  buried  munitions,  the  original  intent  of  this  project  was  to  apply  that  knowledge 
to  the  processing  and  interpretation  of  marine  data  collected  by  other  groups  sponsored  by  SERDP 
and  ESTCP.  However,  the  relevant  hardware  based  projects  have  been  delayed  or  cancelled  and  we 
have  not  yet  obtained  any  data  from  other  projects.  We  expect  that  once  good  quality  marine  data 
are  available  that  additional  research  challenges  will  arise.  At  this  stage  we  believe  that  the  main 
challenges  associated  with  practical  use  of  EMI  sensing  to  the  underwater  environment  are  centered 
around  the  logistics  of  operating  a  sensor  underwater  in  close  proximity  to  the  sea-bottom.  Specific 
challenges  to  address  include: 

1)  Deploying  the  EMI  system  so  that  it  remains  within  about  1  m  of  the  sea-bottom. 

2)  Maintaining  good  detection  performance  at  large  standoffs  in  the  presence  of  conductive  sea¬ 
water.  This  will  require  the  use  of  larger  transmitter  loops  which  in  turn  cause  a  large  conductive 
response  from  the  sea-water.  Either  differential  measurements  or  other  background  subtraction 
methods  would  be  required  for  background  subtraction. 

3)  Accurately  positioning  the  EMI  system  across  a  wide  range  of  water  depths  and  subsurface  con¬ 
ditions  where  existing  methods  based  on  layback  calculations  or  Ultra  Short  Baseline  (USBL) 
acoustics  are  not  always  accurate. 

4)  For  operation  in  a  cued  mode,  there  are  significant  challenges  associated  with  relocating  the 
anomaly  of  interest,  and  physically  positioning  and  holding  the  sensor  over  the  correct  location 
for  the  duration  of  the  measurements. 

As  with  the  terrestrial  environment  a  range  of  different  configurations  of  sensors  will  be  required 
to  address  the  wide  range  of  marine  conditions  that  can  be  encountered. 
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